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ABSTRACT 


This  thesis  is  a  report  of  an  experimental  investigation  of  the 
influence  of  horizontal  mechanically  induced  simple  harmonic  vibrations 
upon  the  rate  of  heat  transfer  by  free  convection  from  a  horizontal  heated 
cylinder  to  air.   The  diameter  of  the  cylinder  was  7/8",  and  the  ranges  of 
the  principal  experimental  variables  were  as  follows: 

temperature  potential,    £,.t:  150°F  i   IF 

frequency  of  vibration,    f :  l6k  cps 

amplitude  of  vibration,    a:  0.025  "to   0.128" 

intensity  of  vibration,    af:  0.21  to   1.11  ft/sec 

The  data  show  that  the   influence  of  vibrations  upon  the  heat- 
transfer  coefficient  is   negligible   for  intensities  of  vibration,    af,    less 
than  0.3  ft/sec;    above   this   so-called  "critical  intensity,"   the  effect  of 
vibrations   is   to   increase  the  heat-transfer  coefficient  significantly. 
The   quantitative  heat  transfer  results  of  this   investigation  show  close 
agreement  with  the  results   reported  for  horizontal  cylinders   in  the  pre- 
sence of  high  intensity  horizontal  sound  fields.     An  elementary  flow- 
visualization  examination  utilizing  smoke  was   performed  which  indicated 
that  the  free   convective  flow  near  the   lower  test  cylinder  stagnation 
point  remained  laminar  in  the  presence  of  vibrations .      The   quantitative  heat 
transfer  results,    supported  by  the   flow  study,    lead  to  the   conclusion  that 
the  physical  mechanism  of  heat  transfer  for  mechanically  induced  horizontal 
vibration  of  a  cylinder  is   the  same  as   for  acoustically  induced  horizontal 
vibrations   in  the   fluid  medium  surrounding  the   cylinder.      The   close   quanti- 
tative  agreement  of  the   results   for  the  two  methods   of  vibration  also   indi- 
cates  that  the   influence  of  the  amplitude  of  vibration  to   cylinder  diameter 
ratio,    a/DQ,   which  differed  by  a  factor  of  12.3,    is  not  significant  for  the 
range  of  a/D0   ratios   considered. 
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CHAPTER  I 
INTRODUCTION 

The   study  of  heat  transfer  in  the   presence  of  vibrations   is   a  rela- 
tively new  area  of  investigation  that  has   received  increasing  attention  in 
recent  years.      Several  experimental  investigations  have   indicated  a  marked 
increase   in  convective  heat-transfer  rates   in  the  presence  of  strong  vibra- 
tions.     Since  heat  transfer  is   a  limiting  factor  in  many  areas  of  engineer- 
ing design,    it   is    important  to  account   for  the  effects  of  vibration  on 
heat  transfer  where  vibration  is  present.      In  certain  cases,   advantage  may 
be   taken  of  increased  heat-transfer  rates  to   render  certain  heat-transfer 
dependent  processes  more  efficient.      As   an  example,   Raben  and  Dietert   ^    ' 
have   conducted  an  investigation  designed  to  use  vibrations   to   increase 
heat-transfer  rates   in  a  saline  water  conversion  process. 

The   influence  of  vibrations  upon  the   free   convective  heat-transfer 
rate   from  a  heated  surface   can  be   investigated  experimentally  by  two  dif- 
ferent methods.      In  one  method  the   heated  surface   is   held  stationary 
within  a  fluid  medium.      An  oscillatory  motion   is   impressed  upon  this   fluid 

medium  --  that  is,  this  method  requires  the  establishment  of  a  sound  field 

1 
in  the   vicinity  of  the   heated  surface   and  will,    therefore,   be   referred  to 

subsequently  as   the    "sound  method."     An  alternate  experimental  approach  is 
to   impress  a  vibratory  motion  on  the  heated  surface   itself,    leaving  the   fluid 
medium  surrounding  the   surface  otherwise   undisturbed.      This   second  method 
generally  employs   a  resonant  mechanical  system  to  achieve   sustained  vibra- 
tions of  the  heat-transfer  surface   and  will,    therefore,   be   referred  to   as 
the    "mechanical  vibrations  method."      The  objective   achieved  by  the   two 
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methods   is   identical.      That   is,    the   relative  velocity  vector  between  the 

fluid  medium  and  the  heated  surface  undergoes   a  periodic  change  of  magnitude 

(2) 
and  sign.      Indeed,    an  analytical  study ^    '   has   shown  that,    for  sufficiently 

long  wave-lengths,    identical  fluid  flows   should  exist  near  a  surface   if  the 

direction,    amplitude,    and  frequency  of  vibration  are  the   same,    regardless  of 

which  of  the   two  methods   is  used  to   create  the  vibration. 

An  early  attempt  to   correlate   the  effects   of  vibration  on  heat  trans- 
fer is   that  of  Martinelli  and  Boelter^'.      In  this   investigation,    a  3 A" 
diameter  horizontal  electrically  heated  tube,    immersed  in  water,   was   sub- 
jected to  vertical  sinusoidal  vibration.     Amplitudes  up  to   0.10"   and  fre- 
quencies up  to     kO  cps  were  employed.      The   range  of  temperature  potential, 
&t,   extended  from  8  to  45°F.     No   influence  on  the  heat-transfer  coefficient 
was  observed  in  the   lower  range  of  vibration  Reynolds  number;    this  result 
was   attributed  to  the   dominance  of  free   convection  in  this   range.     At  higher 
vibration  Reynolds  numbers,    however,    a  gain  of  as  much  as   ^-00  per  cent  rela- 
tive  to   the  non-vibration  case  was  observed.     Unfortunately,   Boelter,    in  a 
communication  to  Lemlich        ,   has   reported  that  these   results   could  not  be 
verified  in  a  later  series   of  experiments . 

The  effect  of  vibrations  on  convective  heat-transfer  rates   from 
horizontal  wires   to   air  has  been  studied  by     Lemlich^    ' .      Three  heated 
wires   of  varying  diameter   (0.0253",    O.O396"   and  0.0810")  were   set  in  vibra- 
tion at  frequencies  between  39  and  122  cps  with  a  maximum  amplitude  of  0.115". 
Increases   in  the  heat-transfer  rates  were   reported  for  increases   in  frequency 
or  amplitude.      Improvements   in  heat-transfer  coefficients  by  as  much  as   a 
factor  of  four  were   reported.      In  addition,   Lemlich  reported  similar  results 
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for  both  vertical  and  horizontal  vibrations   and  concluded  that  the  varia- 
tion in  heat  transfer  was   independent  of  the  direction  of  vibration. 

Tsui^'   has   conducted  an  experimental  investigation  of  the  effect 
of  transverse  vibrations  on  the   heat-transfer  rate   from  a  heated  vertical 
plane   surface.      Increases   in  heat-transfer  coefficients  up  to  25  per  cent 
were  noted.      Tsui  attempted  to  treat  the  system  analytically,   but  conclud- 
ed that  the   general  equations   correlating  vibration  and  heat  transfer  are 
insolvable  with  present  techniques.      Shine ^^   performed  experiments   similar 
to  those  of  Tsui.     A  heated  vertical  plate  was   subjected  to  transverse  vi- 
brations of  11  to   315  cps   at  amplitudes  up  to   0.06l".      Plate  temperature 
range  was   from  131  to  279°F.      Significantly,    these   tests   demonstrated  that 
as  the   amplitude -frequency  product,    af,    (arbitrarily  defined  as  the  vibra- 
tion intensity)  was   increased,    the   heat  transfer  was  not  affected  until  a 
specific  value  of  af  was   reached.      Above   this    "critical"  vibration  inten- 
sity,   af  =  O.O83  ft/sec,    the   rate  of  heat  transfer  increased  markedly. 

(7) 
Freedman      '   has   investigated  the  case  of  a  vertical  heated  cylinder 

vibrating  in  an  axial  direction.      No   change   in  heat-transfer  rate  was 

noted  over  the   range  of  vibration  intensities   achieved   (maximum  vibration 

intensity,    af  =  0.2   ft/sec). 

Fand,   Kaye   and  Teleki^    '    ?'   have   recently  conducted  experimental 

investigations   of  the   influence  of  vertical  mechanical  vibrations 

(Fig.    1-b)   upon  the  heat-transfer  rate  by  free   convection  from  a  heated 

horizontal  cylinder  to  air.   The  diameter  of  the  cylinder  was  7/8" >   and 

the  ranges  of  the  variables  were  as  follows:   temperature  potential,^! t, 

25  to  185  F;  amplitude  of  vibration,  a,  0  to  0.l6";  frequency  of  vibration, 
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f,    5^  to  225  cps;    and  intensity  of  vibration,    af,    0  to   1.22  ft/sec.     These 
authors   report  that  the   sole   controlling  vibrational  variable   is   the  vibra- 
tion intensity,    af.   Further,    the   critical  value  of  af  was   about  0.3  ft /sec. 
A  flow-visualization  study  using  smoke   indicated  that  the   fluid -dynamical 
mechanism  which  causes   the   increase   in  heat  transfer  rate   is   vibrationally 
induced  turbulence . 

An  extensive   investigation  involving  the  sound  method,   that   is,   the 
effect  of  sound  fields  upon  the  heat  transfer  by  free   convection  from  a 
stationary  horizontal  cylinder  to  air,    has  been  reported  by  Fand  and 
Kaye  .      In  their  experimental  arrangement  the   sound  vibration  vector 

was   normal  to  both  the   cylinder  axis   and  the   gravity  field   (Fig.    1-a) . 
The   frequency  range   covered  was   from  1100  to   6120  cps;    the  vibration  in- 
tensity range  was   from  0  to   1.3  ft/sec;    and  the   temperature  potential  was 
varied  from  0  to  250  F.      A  critical  vibration  intensity  of  approximately 
0.3  ft /sec  was  observed;   below  this   intensity  the  heat-transfer  coefficient 
was  unaffected.      The  data  also   showed  that  for  sound  waves  whose  half  wave- 
length was   six  or  more   times   the   test  cylinder  diameter,   the   heat-transfer 
coefficient   is   a  function  of  only  two  variables;    namely,    the   temperature 
potential,    At,    and  the  vibration  intensity,    af .      In  the  tests   above   the 
critical  intensity  the   rate  of  heat  transfer  increased  rapidly  until  an 
af  of  about  0.6  ft/sec  was   reached;    above   this    intensity  level  the   rate  of 
change  of  heat  transfer  was   less   rapid  and  could  be  expressed  by  the   fol- 
lowing algebraic  expression: 


\     =     0.T22 


A  t   (af )2 


5- 
1/3 

(1) 


where  h     is   the   coefficient  of  heat  transfer  in  the  presence  of  vibra- 

v 

tions,  Btu/hr-ft  -°F;  af  is  the  vibration  intensity,  ft /sec;  and  ^\t 
is  the  temperature  potential,   F. 

Fand  and  Kaye  examined  the  physical  mechanism  which  caused  the 

heat-transfer  rate   to   increase   in  the  presence  of  sound  by  means  of  a 
smoke   flow-visualization  technique .      This  study  showed  that  at  the   criti- 
cal sound  intensity  the   typical  free-convective  boundary- layer  flow  pat- 
tern around  the  heated  cylinder  is   disrupted,    and  a  fundamentally  differ- 
ent type  of  boundary-layer  flow,    called  the rmo acoustic  streaming, develops . 
This   thermoacoustic  streaming  appears   as   a  special  type  of  vortex  motion. 
In  the  presence  of  both  horizontal  vibration  and  a  temperature   gradient, 
a  pair  of  vortices   alternately  form  in  the  two   quadrants   above  the  heated 
cylinder  surface.      The   inception  of  this  oscillating  vortex  pair  forma- 
tion occurs   at  the   critical  vibration  intensity,    and  the   further  develop- 
ment of  these  vortices   corresponds   to   that  range  of  af  wherein  the  heat- 
transfer  coefficient  increases  most  rapidly.      A  further  increase   in  the 
vibration  intensity  serves   to  more   fully  develop  these  vortices,   but  does 
not  alter  their  form. 

The   investigations   summarized  above  have,    for  the  most  part,  been 
of  an  experimental  nature.      All  attempts   to  obtain  a  general  analytical 
solution  to   the  problem  have  been  unsuccessful  thus   far.      In  view  of  this, 
it  appears  necessary  to   add  to   the  knowledge  of  the  physical  relationships 
between  vibrations   and  heat  transfer  by  empirical  means .      The  primary 
consideration  involved  in  the  experimental  program  undertaken  in  this 
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work  is   to  establish  the  heat-transfer  "behavior  of  horizontal  cylinders 
under  horizontal  mechanically  induced  sinusoidal  vibration,   where   the 
vibration  vector  and  the   force  of  gravity  are   in  the  same   relationship  to 
a  heated  test  cylinder  as   in  the   sound  investigations   of  Fand  and  Kaye . 
No  direct  comparison  of  the   sound  method  and  the  mechanical  vibrations 
method  with  both  the  test  cylinder  axis   and  the  vibration  vector  normal 
to  the  direction  of  the   gravity  field   (Fig.   1-a)   has  been  reported.     As 
noted  above,    investigations  made  of  the   influence  of  vertical  vibrations 
on  the  heat-transfer  rate  of  a  horizontal  cylinder  to  air   (Fig.    1-b)   have 
indicated  that  the  basic  mechanism  of  convective  air  flow  around  the   cyl- 
inder is  vibrationally  induced  turbulence,   not  thermoacoustic   streaming 
as   observed  in  the   sound  method.      Therefore,    it  was   desired  to  determine 
if  the   character  of  heat  transfer  and  the  associated  thermoacoustic   stream- 
ing phenomenon  exhibited  in  the   sound  method  can  be   related  to  the  mechani- 
cal vibrations  method.      Accordingly,    this   study  is   an  experimental  inves- 
tigation of  the   influence  of  horizontal  mechanically  induced  simple  har- 
monic vibrations  upon  the  heat-transfer  rate  by  free   convection  from  a 
heated  horizontal  cylinder  to   air. 


CHAPTER  II 
GENERAL  DISCUSSION  OF  EXPERIMENTAL  DESIGN 

The  primary  objective  of  this  investigation  was  to  determine  experi- 
mentally the  influence  of  transverse  horizontal  mechanically  induced  simple 
harmonic  vibrations  upon  the  heat-transfer  rate  by  free  convection  from  a 
heated  horizontal  cylinder  to  air.   This  configuration  (Fig.  1-a)  was 
chosen  in  order  to  provide  direct  comparison  with  available  data  for  sta- 
tionary horizontal  cylinders  in  the  presence  of  transverse!  horizontal 
sound  fields.   In  particular,  the  experimental  program  was  designed  to 
establish  whether  or  not  a  quantitative  comparison  of  the  heat-transfer 
variation  for  the  two  cases  at  identical  temperature  potentials,  /it, 
and  intensities  of  vibration,  af,  does  in  fact  exist.   To  this  end,  then, 
the  following  experimental  goals  were  established: 

a.  To  construct  a  mechanical  vibratory  system  which,  at  resonance,  could 
drive  the  test  cylinder  at  vibrational  intensities,  af,  up  to 

1.50  ft/sec 

b.  To  provide  a  suitable  means  to  measure  vibration  amplitude  and 
frequency  accurately 

c.  To  construct  a  heat  transfer  test  cylinder  which  could  operate 
under  the  imposed  vibration  at  surface  temperatures  up  to  350  F 

d.  To  provide  a  power  supply  to  heat  the  test  cylinder  and  a  tem- 
perature measuring  system  for  accurate  temperature  measurement 
of  the  test  cylinder  surface. 
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Based  on  the  experience  of  others^  >  "',  it  was  determined  that  the 
test  cylinder  should  be  rigidly  coupled  to  a  resonant  mechanical  system. 
At  resonance  the  force  magnification  factor  could  be  utilized  for  obtain- 
ing maximum  vibrational  amplitude.  Also,  although  the  test  program  con- 
templated here  involved  only  horizontal  vibrations,  it  was  considered 
economical  to  design  the  vibration  test  equipment  with  a  view  toward 
future  work  employing  other  vibration  directions.   Hence,  it  was  decided 
to  incorporate  into  the  vibration  system  means  by  which  the  vibration 
direction  could  be  varied  from  0  to  90°  with  respect  to  the  direction  of 
the  force  of  gravity. 

The  basic  element  of  the  mechanical  system  chosen  to  provide  vibra- 
tion was   a  simple   free-free  beam  as   shown  in  Fig.   2.      The   free-free  beam 
offers   the   advantage  of  confining  the  major  forces,   which  are  the   inertia 
forces,   to   the  beam  itself  and  not  to.  the   supporting  structure.      Support- 
ing the   free-free  beam  at  the  nodal  points  minimizes   any  hysteresis   losses 
in  the  supports   since   there   is  no   force   transmission  to   the  supports. 
Therefore,    less  energy  absorption  occurs   at  the  nodal  supports   than  would 
occur,    for  example,    at  the   supports   of  a  clamp-clamp  or  hinge -hinge  beam. 

The   test  cylinder  was   rigidly  connected  to  the   test  cylinder  by 
means  of  a  yoke.      The  beam  itself  was   supported  at  the  nodes  by  a  "gimbal" 
arrangement  attached  to   a  rigid  "alignment  platform"    as   shown  in  Figs.   2 
and  3«      The  distance  between  the   gimbal  supports   could  be  varied  to  the 
computed  nodal  positions  of  the  beam  when  the  beam  was   placed  in  resonant 
vibration.      The   alignment  platform  was   attached  to   a  MB  Electronics  Model 
CIF  Vibration  Exciter  which  provided  the  external  driving  force .      The 
motion  of  the  test  cylinder  resulting  from  an  electric  excitation  of  the 
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coupled  dynamic  system  was  a  uniform,  one -dimensional  sinusoidal  trans- 
verse vibration.   The  frequency  of  vibration  employed  was  10^  cps;  however, 
beams  designed  for  other  resonant  frequencies  could  be  used,  as  provision 
was  made  for  fitting  beams  of  variable  dimensions  into  the  vibration  test 
stand.  Vibration  amplitude  was  measured  by  a  cathetometer. 

The  test  cylinder,  which  was  heated  by  an  internal  electrical  heater, 

had  a  diameter  of  7/8" •   Test  cylinders  previously  used  in  other  experi- 

(8,  9) 
ments       involving  heat  transfer  under  conditions  of  vibration  employed 

guard  heaters  on  either  side  of  the  central  test  section  heater  in  order 
to  eliminate  the  problem  of  axial  heat  conduction  away  from  the  test  sec- 
tion.  The  complexity  resulting  from  this  type  of  construction  —  that 
is,  the  use  of  three  interconnecting  parts,  each  with  its  own  heating 
element  and  thermocouple  arrangement  --  has  resulted  in  premature  failure 
by  fatigue.   In  order  to  avoid  fatigue  failure,  the  test  cylinder  for  this 
project  was  designed  with  one  heater;  a  metering  device  was  installed  at 
each  end  of  the  cylinder  to  measure  the  axial  conduction  heat  losses  to 
the  test  cylinder  supports . 

A  total  of  eight  thermocouples  were  installed  in  the  test  cylinder 
for  surface  temperature  measurement.  The  location  of  these  thermocouples 
is  as  shown  in  Fig.  5*  Each  copper-constantan  thermocouple  was  encapsu- 
lated in  the  end  of  a  copper-plated  stainless  steel  tube.  The  tube  was 
then  closely  fitted  into  its  individual  hole  drilled  axially  into  the 
cylinder  wall.   Thus  the  thermocouples  were  firmly  fixed  in  position, 
and  a  good  thermal  contact  between  the  junctions  and  the  cylinder  was 
established. 
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There  vas   an  obvious   danger  of  test  cylinder  lead-wire  "breakage 
under  vib rationally  induced  loads   as   high  as   200  g's   during  the  course 
of  the   tests.      Several  special  techniques  were  used  to   avoid  wire  break- 
age  and  are   described  in  detail  in  Chapter  V,   Section  1-c . 

The  vibration  test  stand  itself  was  located  in  a  special  test  room 
8  x  13  x  15 '  high.   A  k   x  6'  glass  window  was  placed  in  the  wall  between 
the  test  stand  and  the  adjacent  control  room.   Thus,  optical  amplitude 
measurements  and  all  other  test  operations  could  be  conducted  without 
disturbing  the  air  in  the  vicinity  of  the  test  cylinder. 
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CHAPTER  III 
APPARATUS 

The  apparatus  can  "be  conveniently  divided  into  two  basic  categories 
with  subdivisions  as  follows: 

1.  Vibration  Apparatus 

1-a.  General 

1-b .  Resonant  beam 

1-c .  Test  stand 

1-d.  Vibration  exciter  and  exciter  power  supply 

1-e .  Amplitude  and  frequency  measuring  apparatus 

2.  Heat  Transfer  Apparatus 

2-a.  General 

2-b.  Test  cylinder 

2-c .  Power  supply 

2-d.  Temperature  measuring  apparatus 

1.  Vibration  Apparatus 

1-a.   General 

In  order  to  produce  the  intensities  of  vibration  required  for  this 
investigation,  it  was  necessary  to  employ  a  mechanical  resonant  system. 
This  was  evident  from  the  characteristics  of  the  available  vibration 
exciter;  the  maximum  acceleration  obtainable  from  this  exciter,  using 
direct  force  take-off,  was  about  75  g's.  As  stated  in  Chapter  II,  ac- 
celeration loads  up  to  approximately  200  g's  were  required  to  obtain  the 
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desired  value  of  af  =  1.5.   To  obtain  these  high  accelerations,  it  was 
necessary  to  interpose  a  mechanical  resonant  system  between  the  vibration 
exciter  and  the  test  cylinder.   This  is  the  general  technique  described 
in  the  literature  used  to  produce  high  values  of  vibration  intensity. 

An  additional  requirement  of  the  vibration  system  was  that  the  vi- 
bration motion  produced  had  to  be  sinusoidal  and  one -dimensional  in 
character.   This  requirement  restricted  the  selection  of  the  resonant 
system  to  one  which  would  vibrate  essentially  in  its  first  natural  mode 
with  negligible  effect  from  any  higher  harmonics.  A  prismatical  beam  was 
the  system  chosen  which  satisfied  both  the  resonance  and  sinusoidal,  one- 
dimensional  motion  requirements. 

Other  mechanical  resonant  systems  could  have  been  devised  to  pro- 
duce the  motion  required  --  such  as  a  system  using  heavy  coiled  springs, 
for  example.   However,  the  simple  geometry  of  the  rectangular  beam  lends 
itself  more  conveniently  to  accurate  analysis,  design  and  construction. 

Resonant  beam  vibration  with  clamp-clamp,  hinge-hinge  and  free -free 
boundary  conditions  were  investigated.   The  free-free  vibration  of  a  beam 
in  the  first  mode  was  selected  as  the  preferable  method  for  obtaining  the 
required  vibration  intensity  because  this  system  exhibits  the  least  damp- 
ing.  In  free-free  beam  vibration  the  major  forces,  which  are  the  inertia 
forces,  are  balanced  within  the  beam  itself  and  not  transmitted  to  the 
structure  supporting  the  beam.   Thus,  the  joint  type  of  damping  loss  is 
eliminated.  An  additional  advantage  results  from  this  absence  of  major 
force  transmission  at  the  supports:   the  supporting  structure  is  largely 
free  of  induced  vibrations  and  can  be  designed  accordingly. 
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One  disadvantage   in  employing  a  free -free   resonant  beam  compared  to 
a  clamp-clamp  resonant  beam  is   the  necessity  to  provide   a  different  beam 
for  each  vibration  frequency  desired.      However,   since  it  has  been  demon- 
strated by  other  investigators^    '    *'   that  the   sole  vibrational  variable 
affecting  heat-transfer  rates   is   the   intensity  of  vibration  —  the  ampli- 
tude-frequency product,    af.      Hence,    it  is   possible   to   investigate  the 
influence  of  vibration  intensity  on  heat  transfer  by  varying  the   amplitude 
of  vibration  at  constant  frequency. 

1-b .      Resonant  beam 

The   complete   analysis  of  the   free-free  beam  with  center-span  test 
cylinder  load  is   given  in  detail  in  Appendix  A.      The  results   and  interpre- 
tation of  this   analysis   are  presented  here. 

In  the  analysis   it  was   assumed  that  the  beam  at  its   resonant  fre- 
quency,  has   a  motion  essentially  of  the  first  mode  only,    the  effects   of 
the  higher  harmonics  being  negligible.      This   assumption  is   considered 
valid  for  a  uniform,    carefully  machined,   prismatical  beam.      Further,    tests 
showed  that  the  vibratory  motion  of  the  beam  was,    in  fact,   sinusoidal  and 
one -dimensional  in  character. 

The  resonant  frequency  equation  of  the  system,    as   derived  in  Appen- 
dix A,    in  terms  of  beam  mass,   m^,    and  the  mass   at  beam  center-span,   mc, 
is: 

™b  Q(a     (sechof/  +   cosoCa) 

n£  ( s inflf/  +   cos c{J(  tan  h tf/  )  (2) 

where  A     is  one-half  the  beam  length  and  0(    is  defined  as: 

k  _  i2_r  m2 

^        "     Eg.        W 
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p 

E   =   modulus  of  elasticity,  lb/in 
g   =   gravitational  constant,  in/sec 
X  =   "beam  weight  density,  lb/in 
00  =   circular  frequency,  rad/sec 
h   =   beam  height  in  motion  direction,  in 
Equation  (2)  is  plotted  in  Fig.  6. 

In  terms  of  the  beam  physical  parameters,  the  resonant  frequency 
can  be  expressed  as: 

1      &  1/2    E   1/2    0 

Therefore,  for  a  beam  of  given  dimensions,  the  vibration  intensity  is: 

E  2/2 
af  =  Constant  x  a  x  (— )  (k) 

Consideration  of  maximum  bending  stress  leads  to  the  following  rela- 
tionship between  maximum  bending  stress  in  the  beam,  (T_. ov,    and  amplitude 
of  vibration,  a,  as  derived  in  Appendix  A  and  plotted  in  Fig.  9*. 

^"max  E       h        ,  „  *,  0  \  r *-\ 

T~     =     2       T2   ^a  function  of  ^  )  (5) 

Equations    (2),    (4)   and   (5)   form  the  essential  basis   for  resonant 
free -free  beam  selection.     Fig.    9  clearly  shows   that  a  high  value  of  &y   is 
desired  in  order  to  minimize  bending  stress   in  the  beam,    and  Fig.    6  indicates 
that  large  values  of    Q(J(     correspond  to  large  values  of  the  beam  mass   to 
center  load  mass   ratio,   m.  /m   .      The  weight  of  the  center  load  was   known  to 
be   approximately  5-5  lbs.      This  weight  included  the  weights  of  all  the 
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moving  parts  attached  to  the  center  of  the  "beam:   the  test  cylinder  yoke, 
the  moving  element  of  the  vibration  exciter,  and  the  parts  'used  to  couple 
the  vibrator  moving  element  and  test  cylinder  yoke  to  the  beam.  There- 
fore, it  is  clear  from  Fig.  6  that  for  a  reasonable  beam  size,  the  value 
of  (Xa       would  have  to  be  within  the  range  of  about  2.20  to  2.28.  These 
values  of  Q(j[    represent  beam  weights  from  19  lbs  to  52  lbs  as  obtained 
from  Fig.  6,  using  a  value  for  center-span  load  of  5*5  lbs.   In  the  actual 
beam  selection  a  compromise  value  of  0(4  -     2.25  and  a.  corresponding  beam 
weight  of  31*6  pounds  were  used  for  purposes  of  analysis.  A  lower  value 
of  (Xa     would  impose  a  greater  bending  stress,  and  a  higher  value  of  0(y{ 
would  require  a  beam  larger  than  could  be  conveniently  used. 

Equation  (k)   indicates  that  a  high  value  of  E/y   is  desired. 

Additional  considerations  involved  in  the  beam  design  include  the> 
vibration  frequency  and  fatigue  strength.   The  frequency  at  which  reson- 
ance was  desired  was  in  the  neighborhood  of  100  cps .   This  frequency  was 
selected  as  being  high  enough  to  be  well  within  the  frequency  ranges  used 
by  other  investigators  to  demonstrate  marked  effects  of  vibration  on  heat- 
transfer  rates,  and  low  enough  to  avoid  more  severe  inertia  loads  at  yet 
higher  frequencies.   The  loading  is  proportional  to  the  square  of  the 
frequency,  and  even  at  frequencies  of  100  cps,  inertia  loads  of  about 
200  g's  were  anticipated.   The  vibration  test  stand  was  designed  so  that 
beams  of  varying  dimensions  could  be  installed;  hence,  frequency  could  be 
varied  by  changing  beam  size. 

Fatigue  strength  was  a  primary  factor  involved  in  beam  material 
selection.   The  beam  was  required  to  sustain  complete  stress  reversals 
for  long  periods  of  time .  Beam  failure  might  have  resulted  in  test 
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cylinder  destruction  or,    at  least,   would  have   caused  considerable  delay 
in  the   test  program. 

Of  the   several  beam  materials   considered,    aluminum  and  steel  were 
given  consideration  based  on  strength  and  fatigue   characteristics,    avail- 
ability,  machinability,    and  the   requirement  of  Equation   (k) ,   namely,    a 
high  E/v<    ratio.      From  the   stress   standpoint,    the  best  available   aluminum 
alloy  was   202^-T^,    a  special  heat-treated  alloy  which  has  a  yield  stress 

of  50*000  lb/in   ,    and  the  best  available   steel  alloy  was   a  special  treat- 

o 
ment  steel,    "HY-100,"  which  has   a  yield  stress   in  excess  of  100,000  lb/in    . 

The  endurance   limit   for  Al  2024-T^  is   quoted  as   20,000  lb/in2.      This 
is   an  arbitrarily  defined  endurance   limit,    for  it   is  well  known  that  the 
fatigue   characteristics   of  non-ferrous  metals   do  not  exhibit  a  limiting 
value  of  fatigue   strength  below  which  the  material  strength  is   independ- 
ent of  the  number  of  stress   cycles.     No   fatigue   data  for  HY-100  steel  is 
currently  available.      It   is   known  that  the  endurance   limit   for  most 
alloy  steels    is   approximately  one-half  the   yield  strength;    therefore, 
"HY-100"   steel  was   considered  to  be   adequate   from  the   fatigue   standpoint 
up  to   stress   levels   of  approximately  50,000  lb/in    . 

Each  of  these   two  metals   could  be  machined  readily;    and  aluminum 
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and  steel  have   about  the   same  E/^      ratios    --  1.0  x  10     and  1.1  x  10 

inches,    respectively  --  giving  each  material  about  the   same   footing  rela- 
tive  to  the   requirement  of  Equation   (k) . 

The  width,   b,    of  the  beam  was   established  as   1.50"   initially  in  order 
that  existing  test   cylinder  support  equipment   could  be   used.      In  a  study 
made  of  Al  2024-TU  and  HY-100  beams,    the   first  comparison  was   made  on  the 
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"basis   of  equivalent  "beam  size,    and  the   second  comparison  was  made   for 
identical  values  of    (Xa      .      The   results   of  the   equal  beam  size   calcula- 
tions made   in  accordance  with  the  method  given  in  Appendix  B  were: 

Beam  Size    --   I.5O11   X  l.kk"   x   52,00" 

<A$  f(cps)  af   (ft/sec)      (T^  (lb/in2) 

Steel  2.25  100  1.50  52,000 

Aluminum  2.15  100  I.50  25,500 

The   results  of  the   comparison  on  the  basis   of  identical  values   of 
were: 

Size                          of/  f(cps)        af   (ft/sec)    (T^  (lb/in2) 

Steel            1.50"xl.Wx52.00"        2.25  100                   I.50                     52,000 

Aluminum     1.50"x2.90"x73.00"        2.25  100                   1.50                     18,000 

These   figures   indicate,    for    (Xa    =     2.25,    that  either  material  would 
be  satisfactory.      Considering  the  superior  fatigue   characteristics   of 
steel  as   compared  to   aluminum,    however,   HY-100  was   selected  as   the  beam 
material  for  the   final  design. 

The  HY-100  obtained  for  the  project  had  the  following  alloy  content 
(MIL-S-2015UA): 

Carbon  Nickel  Chromium  Molybdenum 

0.20-0.21  3.5O  1.50  0.1+0 
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Two   standard  tensile   specimen  tests   gave   the   following  results: 

No.    1  No.    2 

Yield  strength   (lb/in2)  106,000  102:, 200 

Ultimate  strength   (lb/in2)  124,300  111,900 

The  dimensions   and  characteristics  of  the  "beam  used  in  the  test 
stand  were   as   follows: 

width,  b  height,    h  length,    21  f   (calc.)  f   (actual)  (X 


J 


1.47V'  1A0811  50.00"  105.3  cps  104cps  2.2^3 

Small  variations   from  the  planned  dimensions  were  necessary  because  of 
the  machining  required  to  produce   a  straight  beam  with  a  63  micro-inch 
finish  from  the   available   stock.      In  addition,    the   length  of  the  beam  was 
adjusted  to   raise   the   calculated  resonant  frequency  somewhat  above   100  cps 
to   allow     for  the   added  mass   of  the   gimbals  mounted  at  the  beam  nodes   and 
because   the  weight  placed  at  beam  center  span  was   not  accurately  known 
until  all  parts  were   constructed.      The  measured  resonant  frequency  was 
about  1  per  cent  less   than  the   calculated  value. 

1-C.      Test  stand 

The  vibration  test  stand  and  its  components  are  shown  in  Figs.  2, 
3  and  k.     The  purpose  of  the  test  stand  was  to  provide  a  means  of  support 
for  the  resonant  beam  and  a  method  for  accurate  alignment  of  the  resonant 
beam  with  respect  to  the  vibration  exciter. 

Resonant  beam  support  at  the  nodes  of  vibration  was  accomplished  by 
two  gimbal  assemblies.  At  each  node  an  aluminum  gimbal  frame  was  fitted 
to  the  resonant  beam  by  means  of  a  stainless  steel  slipper.   The  beam 
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slipper  was  incorporated  into  the  gimbal  design  so  that  beams  of  different 
sizes  could  be  fitted  into  the  gimbal  frame  --  that  is,  each  beam  with  its 
own  pair  of  slippers  could  rest  in  the  existing  gimbal  frames.  A  beam 
clamp  plate  served  to  hold  the  beam  firmly  within  the  gimbal  frame  by  means 
of  two  l/k"   set  screws  and  lock-nuts.   To  allow  for  gimbal  frame  oscilla- 
tory rotation  caused  by  beam  deflection  during  vibration,  hardened  steel 
journals  were  fitted  into  both  sides  of  the  gimbal  frame  at  its  axis  of 
rotation.   The  journals  were  first  put  through  the  holes  in  the  gimbal 
supports  and  then  seated  in  the  gimbal  frame  holes  with  5/l6"  shoulder 
bolts .   Operation  of  the  apparatus  later  indicated  that  the  journals  came 
loose  in  the  gimbal  frame  holes,  causing  local  heating  and  eventual  seiz- 
ure of  the  frame  within  the  gimbal  supports.   To  correct  this  situation, 
each  journal  was  keyed  to  the  gimbal  frame,  effectively  preventing  any 
relative  motion  between  the  journal  and  gimbal  frame .   In  addition,  a 
journal  oil  supply  system  was  installed;  oil  cups  were  mounted  on  top  of 
each  gimbal  support  post  with  oil  wicks  leading  down  to  the  journal  through 
holes  in  these  support  posts.   To  distribute  the  oil  around  the  journal,  a 
circumferential  l/l6"  x  l/l6"  groove  was  milled  in  the  center  of  the  9/l6" 
diameter  hole  of  each  gimbal  support  post.  A  10  per  cent  colloidal  dis- 
persion of  molybdenum-disulfide  in  oil  was  used  as  the  lubricant. 

Each  gimbal  assembly  was  secured  to  an  alignment  platform  by  four 
3/8"  bolts.   Since  the  nodal  position  was  located  by  calculation  (see 
Appendix  A),  small  errors  in  gimbal  frame  mounting  position  did  exist. 
Thus,  some  vertical  force  was  applied  to  the  gimbal  supports.  A  rigid 
connection  of  the  gimbal  assembly  to  the  alignment  platform  would  have 
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introduced  a  significant  amount  of  damping.   Therefore,  provision  for 
small  relative  motion  between  the  gimbal  assembly  and  the  alignment  plat- 
form was  made  by  placing  a  l/32"  sheet  of  neoprene  rubber  between  them. 

The  alignment  platform  was  designed  to  insure  accurate  alignment  of 
the  resonant  beam  with  respect  to  the  vibration  exciter  driver  axis.  Ac- 
curate positioning  of  the  resonant  beam  was  necessary  so  that  no  lateral 
force  component  could  be  applied  to  the  exciter  driver  coil  through  the 
exciter  extension.   Clearance  between  the  exciter  driver  and  field  coils 
was  such  that  any  misalignment  of  the  exciter  extension  could  cause  the 
driver  coil  to  come  in  damaging  contact  with  the  field  coil.   The  align- 
ment platform,  machined  from  a  standard  8"  structural  channel  section, 
was  slotted  along  the  top  center-line  to  receive  the  pad  projections  of 
the  gimbal  supports .   This  arrangement  prevented  any  lateral  motion  of  the 
resonant  beam.   To  facilitate  moving  the  gimbal  assemblies  to  the  resonant 
beam  node  positions  four  slots,  16"  in  length,  were  machined  in  the  align- 
ment platform. 

Small  stainless  steel  bench  marks  were  located  on  the  alignment 
platform  on  either  side  of  the  central  1  3/8"  diameter  hole  for  the  alumi- 
num exciter  extension.   These  bench  marks  served  as  a  convenient  aid  in 
locating  the  gimbal  assemblies  at  equal  distances  from  the  beam  center. 

A  special  brass  alignment  plug  was  made  to  fit  into  the  1  3/°" 
diameter  hole  through  which  the  exciter  extension  was  centered.   The 
alignment  plug  had  its  center  drilled  out  to  the  exact  diameter  of  the 
exciter  extension.   Thus,  the  installation  position  of  the  alignment 
platform  on  the  vibration  exciter  was  accurately  determined  when  the 
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alignment  plug,   with  the  exciter  extension  through  its   center,   would  just 
slide   into  the   1  3/8"   diameter  hole   in  the   center  of  the   alignment  platform. 

The   alignment  platform  was   attached  to   the   vibration  exciter  by- 
quadrants   and  brackets   in  an  arrangement  which  permitted  rotation  of  the 
entire  exciter-alignment  platform  portion  of  the   test  stand  relative  to 
the  exciter  base.      In  this  manner  the  direction  of  vibration  could  be 
varied  0  to   90     relative   to   the   direction  of  gravity.      The  brackets  were 
bolted  to   the  top  face  of  the  exciter  casting,    and  the   quadrants  were 
fitted  to   the   same  trunnions  which  normally  connect  the  exciter  to   its 
base. 

The  exciter  base  was   in  turn  bolted  to   a  2'    section  of  Ik"  x  Ik  1/2" 
wide-flange  beam  which  served  as   the  base  support  for  the  entire  vibration 
test  stand,    as   shown  in  the  photograph,   Fig.    10.      Adjustable   legs  were 
installed  at  each  corner  of  the  wide-flange  base  to  provide   a  convenient 
means   of  leveling  the  test  stand.      These   legs   supported  the  base  through 
k"   square  pads  of  heavy  rubber  which  prevented  vibration  transmission  to 
the  test  room  floor. 

The  yoke  which  held  the  test  cylinder  at  mid-span  of  the   resonant 
beam  was   a  U-shaped,    finned  structure  which  combined  high  rigidity  with 
minimum  weight.      The  base  plate  of  the  test  cylinder  yoke  was   held  to  the 
resonant  beam  by  means  of  a  steel  clamping  plate.      The  resonant  beam  was 
firmly  clamped  between  these  two  plates   held  together  by  four  l/k"   steel 
bolts.      A  calculation  of  stresses   showed  that  the  maximum  stresses   in  these 
bolts   during  vibration  were   at  all  times  well  below  the   allowable  limit. 
To  prevent   any  corrosive  or  cutting  action  between  the  plates   and  the 
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resonant  beam,    thin  pieces  of  Teflon  sheet  were  placed  "between  the  plates 
and  the  beam. 

The   lead  wires  of  the  test  cylinder  were   fastened  to   the  test  cyl- 
inder in  the  manner  described  in  section  2-b  of  this   chapter.      These 
wires  were  brought  down  along  the  test  cylinder  yoke   and  taped  to  the 
resonant  beam  as   seen   in  Fig.    11.      The   lead  wires  were   held  firmly  along 
the  yoke  by  taping  them  to   a  stiff  U-shaped  wire  which  was   anchored  to 
the  yoke   at  both  ends   and  at  the   center  of  the   resonant  beam  as   can  be 
seen  in  the   close-up;  photograph,   Fig.    11.      The  wires  were   then  taped 
along  the   side  of  the  beam  and  led  off  the  beam  near  the  nodes   to   cable 
clamps   fitted  with  rubber  grommets   on  the   girribal  support  posts .      This 
arrangement  prevented  any  large  relative  motion  between  adjacent  parts 
of  the  wires,    and  lead  wire  breakage  was   avoided  completely. 

1-d.  Vibration  exciter  and  exciter  power  supply 

The  external  driving  force  for  vibration  at  resonance  was  supplied 
by  a  MB  Electronics  Model  C1F  electro dynamic  vibration  exciter.   The 
sinusoidal  force  output  of  this  exciter  was  50  lbs  maximum  within  a  freq- 
uency range  of  6  to  350  cps .   Driver  coil  motion  was  limited  to  a  double 
amplitude  of  l/2";  an  over-travel  limit  switch  was  provided  to  open  the 
driver  coil  circuit  in  case  this  limit  was  exceeded.   The  force  generated 
was  directly  proportional  to  the  value  of  alternating  current  in  the 
driver  coil  for  constant  direct  current  input  to  the  field  coil.  The 
force  rating  of  the  exciter  with  1.75  amperes  input  to  the  field  coil 
was  10  lbs  vector  force  for  each  1.0  ampere  r.m.s.  of  alternating  current 
applied  to  the  driver  coil.   The  maximum  driver  coil  current  was  5-0 
amperes . 
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The  power  supply  for  the  vibration  exciter  was  the  associated  MB 
Electronics  Model  T112031  Power  Supply.  This  power  console  contained  an 
audio-oscillator  with  amplifier  which  supplied  a  sinusoidal  wave  form 
current  to  the  driver  coil  and  included  built-in  impedance  matching  con- 
trols .   The  field  coil  supply  in  the  console  was  a  rectifier  which  pro- 
vided direct  current  up  to  the  rated  1.75  amperes  for  the  exciter  field 
coil. 

l-»e.  Amplitude  and  frequency  measuring  apparatus 

The  vibration  amplitude  of  the  test  cylinder  was  measured  optically 
by  sighting  through  a  cathetometer  on  a  specially  made  target.  The  tar- 
get consisted  of  a  short  piece  of  0.001"  diameter  wire  mounted  on  a 
black  background.  When  the  target  was  illuminated  by  an  intense  spotlight, 
the  image  visible  through  the  cathetometer  microscope  was  a  thin  bright 
line  corresponding  to  the  sharply  illuminated  crest  of  the  wire.  When 
the  cylinder  was  put  into  vibration,  this  image  became  a  bright  distinct 
band  with  a  width  equal  to  the  double  amplitude  of  vibration.  The  ac- 
curacy of  amplitude  measurement  was  better  than  0.002". 

Frequency  measurement  was  obtained  by  feeding  a  signal  from  the 
audio-oscillator  in  the  exciter  driver  coil  circuit  into  a  Hewlett  Packard 
Model  521A  Electronic  Counter.  The  counter  frequency  reading  was  always 
in  good  agreement  with  the  calibrated  dial  of  the  audio -oscillator;  the 
accuracy  of  the  frequency  measurement  was  within  ±   1  per  cent. 
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2.      Heat  Transfer  Apparatus 

2-a.   General 

The  heat  transfer  apparatus  consisted  of  a  test  cylinder,  a  power 
supply  to  heat  the  test  cylinder  and  a  system  of  thermocouples  installed 
in  the  test  cylinder  for  surface  temperature  measurement.   Special  pro- 
blems were  involved  in  the  design  and  construction  of  this  apparatus  due 
to  the  high  inertia  loads  and  temperatures  imposed  on  the  test  cylinder 
during  long  periods  of  vibration. 

Test  cylinders  previously  used  in  other  experiments  involving  heat 
transfer  from  vibrating  cylinders  have  employed  guard  heaters  on  either 
side  of  the  central  test  section  heater  in  order  to  prevent  axial  heat 
conduction  away  from  the  test  section.   The  complexity  of  this  arrange- 
ment —  that  is,  the  use  of  three  interconnecting  sections,  each  with  its 
own  heating  element  and  thermocouples  —  has  resulted  in  premature  fail- 
ure by  fatigue . 

In  order  to  avoid  fatigue  failure,  it  was  decided  to  reduce  the 
complexity  of  construction  of  the  test  cylinder.   This  was  accomplished 
by  placing  only  one  heater  within  the  test  cylinder.   Since  guard  heaters 
were  not  used  to  prevent  axial  heat  conduction  away  from  the  test  sec- 
tion, a  metering  device  was  installed  at  each  end  of  the  cylinder  to  ac- 
curately measure  the  axial  conduction  heat  losses  to  the  test  cylinder 
supports . 

The  test  cylinder  thermocouple  system  was  designed  to  permit  re- 
placement of  the  individual  thermocouples  should  a  failure  occur.  A 
detailed  description  of  the  heat  cylinder  and  its  components  is  given 
below. 
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2-b.   Test  cylinder 

The  test  cylinder,  shown  in  detail  in  Fig.  5>  consisted  of  a  central 
electrically  heated  test  section  which  extended  between  end  supports  that 
were  designed  to  fit  in  the  test  cylinder  yoke  clamps .   The  test  cylinder, 
7/8"  in  diameter,  had  an  overall  length  of  8  l/2"  and  weighed  1.06  lbs. 
The  central  section,  called  the  test  section,  was  5 '20"  long. 

The  electric  heater,    installed  within  the  test  section,   was   made  by 
winding  gage  No.    32  Nichrome   V  wire    (nominal  resistance,    10.26  ohms  per 
foot)  within  a  helical  groove  machined  in  the  surface  of  a  l/2"   diameter 
ceramic  base  .      The  helical  groove  had  a  pitch  of  20  turns  per  inch  and 
served  to   insulate   the  heater  wire.      The  heater  leads   consisted  of  two 
pieces  of  O.O85"   diameter  drill  rod  which  were   centered  at  the  ends   of 
the  ceramic  base.      The  heater  wires  were   attached  to   these  drill  rod  leads 
with  a  special  medium  temperature   solder.      Regular  half-and-half  solder 
could  not  be  used  because  of  the  high  temperatures   involved;    silver  solder 
was   impractical  because  brazing  temperatures    in  excess   of  1000  F  would 
burn  the  No.    32  Nichrome  V  heater  wire.      The   solder  used  was   "TEC"   solder, 
manufactured  by  the   Handy  &  Harman  Brazing  Products   Company.      This   solder, 
composed  of   5  per  cent  silver  and  95  per  cent  cadmium,    has   a  melting 
point  of  6k0  F   --   about  200°F  above  that  of  the   available   lead-tin  solders. 

The  heater  tube  was   closely  fitted  inside   a  copper  sleeve,   O.D.    5/8", 
and  held  in  its   axial  position  by  Teflon  plugs   through  which  the  drill 
rod  heater  leads   extended  as   is   clearly  shown  in  Fig.    5*      This   assembly, 
with  an  O.D.   of  5/8",   was   then  fitted  into   an  aluminum  sleeve  with  an 
O.D.   of  7/8"  which  was   the  outside  element  of  the   test  section.      The 
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copper  sleeve  between  the  ceramic  heater  base  and  the  aluminum  walls  of 
the  test  section  provided  for  a  more  uniform  axial  temperature  distribu- 
tion along  the  test  section. 

Section  X-X  of  Fig.    5  shows   the  detail  of  the   test  cylinder  end 
supports   and  the   copper-Teflon-copper  concentric   rings  which  formed  the 
heat  meter  for  measuring  heat   loss  by  conduction  to  the  test  cylinder 
yoke   supports .      These   concentric   rings   served  as   a  heat  meter  in  the 
following  manner.      The   inner  copper  ring  was   at  the  temperature  of  the 
test  cylinder  end;    similarly,    the  outer  copper  ring  was   at  the  temperature 
of  the   test  cylinder  yoke   support.      The   temperature  distribution  along  the 
copper  rings,    a  length  of  1  3/8",   was   assumed  to  be   constant  because  of 
the  excellent  thermal  conductivity  of  copper.      The  Teflon  insulating  ring 
produced  a  temperature  drop  between  the  two   copper  rings .      This   tempera- 
ture differential  was  measured  by  a  differential  thermocouple  made  by 
soldering  a  gage  No .    2k  constantan  wire  between  the   two   copper  rings   as 
shown  in  Section  X-X,   Fig.    5«      The  heat  loss   from  the   test  cylinder  by 
conduction  to   the  yoke   supports   necessarily  passed  through  these  heat 
meters.      This   heat  loss,   proportional  to  the  temperature   differential, 
^t^im,   was  measured  in  the  manner  described  in  Chapter  V. 

The  copper  and  Teflon  rings  were  assembled  on  the  test  cylinder  ends 
by  press-fitting  them  firmly  in  place.  Each  heat  meter  thus  formed  a  sup- 
port at  the  test  cylinder  end  which  was  held  in  place  by  the  test  cylinder 
yoke   clamp.      Fig.    11  provides   a  clear  view  of  the   final  assembly. 

The  ends  of  the  test  cylinder  were  sealed  by  Teflon  terminal  plates 
which  were  fitted  into  each  end  of  the  aluminum  sleeve  a  distance  of 


27- 

1  l/8" .   Small  Teflon  washers,  l/l6"  thick,  separated  the  Teflon  terminal 
plates  and  the  Teflon  plugs  in  the  ends  of  the  copper  sleeve  surrounding 
the  heater  tube .   The  outer  2  l/l6"  diameter  flange  of  each  terminal 
plate,  shown  in  Fig.  11,  was  the  surface  on  which  all  electrical  lead 
connections  were  made.   The  heater  lead  drill  rods  were  supported  by  l/2" 
long  stainless  steel  plugs  which  were  threaded  into  each  Teflon  terminal 
plate.   The  heater  lead  connections  from  the  power  supply  were  secured  to 
these  plugs  by  means  of  special  nuts  through  which  set  screws  were  placed 
so  as  to  bear  against  the  drill  rod  heater  leads,  thus  insuring  good 
electrical  contact. 

The  measurement  of  test  cylinder  surface  temperature  was  accomplished 
by  means  of  thermocouples  placed  inside  the  aluminum  sleeve  walls.   The 
test  section  was  provided  with  six  thermocouples  arranged  as  indicated  in 
Fig.  5  a"t  Sections  B-B  through  F-F.   Two  additional  thermocouples  were 
installed  under  the  inner  copper  ring  of  each  sandwich  at  Sections  A-A 
and  G-G.   Thus,  these  two  thermocouples  measured  the  temperature  of  the 
inside  copper  ring  within  a  close  approximation. 

The   thermocouples  were  made  of  gage  No.    30  dual  conductor  copper 
constantan  wire  encased  in  fiberglass   insulation.     Each  thermocouple   junc- 
tion was  made  by  removing  the   insulation  from  both  wires   along  about  1" 
at  the  mid-span  of  a  piece  of  thermocouple  wire   approximately  12"   long. 
The   two  wires   at  the   clean  section  were  then  twisted  tightly  together. 
Next,   one  end  of  the   insulated  portion  of  the  wire  was   carefully  pulled 
through  a  section  of  O.O65"   diameter  stainless   steel  tubing  which  had  been 
copper-plated  on  its  outside   surface  to   a  diameter  of  0.068".     When  the 
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twisted  junction  was  pulled  just  inside  this  tube,  Handy  &  Harman  TEC 
solder  was  used  to  encapsulate  the  junction  at  the  tip  of  the  tube.  Each 
tube  was  then  carefully  fitted  into  its  individual  hole  in  the  aluminum 
sleeve  wall.   These  holes,  0.068"  in  diameter,  received  the  thermocouple 
tubes  in  a  tight  fit;  the  soft  copper-plating  on  the  surface  of  the  tubes 
provided  excellent  thermal  contact  between  the  tubes  and  the  test  cylin- 
der wall. 

The  lead  wires  of  the  thermocouples  and  the  two  copper  wires  used 
for  each  heat  meter  described  above  were  closely  fitted  through  holes 
drilled  in  each  Teflon  terminal  plate.  From  these  holes  the  wires  were 
led  to  individual  terminal  affixed  to  the  terminal  plate  by  No.  3  stain- 
less steel  screws  which  were  threaded  into  the  Teflon.  Each  wire  end 
was  formed  into  a  spiral  similar  in  form  to  the  hair-spring  of  a  watch. 
The  spirals  were  secured  with  small  copper  strips  held  by  the  No.  3 
screws.   The  spiral  end  prevented  any  wire  from  crossing  itself  and  break- 
ing from  the  contact  when  under  intense  vibration.   The  lead  wires  lead- 
ing down  the  test  cylinder  yoke  (described  in  Section  1-c  of  this  chapter) 
were  attached  in  a  similar  manner.   Thus,  each  wire  terminal  consisted 
of  spiral  end  from  the  thermocouple  junction,  a  copper  strip,  the  spiral 
end  from  the  lead  wire,  a  second  copper  strip  and,  finally,  two  No.  3 
screws  holding  the  terminal  to  the  Teflon  terminal  plate.   To  prevent 
these  screws  from  coming  loose  during  vibration,  their  slots  were  aligned 
and  small  pieces  of  wire  were  soldered  across  the  slots,  preventing  ro- 
tation of  the  screws. 
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The  extreme  care  used  in  the  installation  and  wiring  of  the  thermo- 
couples proved  worthwhile;  not  one  thermocouple  failure  occurred  during 
the  test  program.  However,  an  attractive  feature  of  this  thermocouple 
installation  method  was  that  an  individual  thermocouple  could  "be  removed 
and  replaced  in  the  event  of  failure . 

The  surface  of  the  test  section  was  polished  to  a  mirror  finish  to 
insure  a  low  thermal  emissivity.  A  low  value  of  emissivity  made  the 
radiation  heat  loss  from  the  test  section  small  compared  to  the  heat 
transferred  by  convection.   Consequently,  errors  made  in  the  calculation 
of  radiation  losses  did  not  introduce  appreciable  errors  in  the  calcula- 
tion of  the  convective  heat  transfer  coefficient.   The  emissivity,  as 
determined  by  comparison  with  a  similar  heat  transfer  test  cylinder  of 
known  emissivity,  was  e  =  0.0^5. 

2-c.   Power  supply 

The  power  supply  to  the  test  cylinder  heater  was  drawn  from  a  500 
watt  constant  voltage  regulator.   Two  variacs,  placed  in  series  with  the 
voltage  regulator  output,  served  for  the  coarse  adjustment  of  the  voltage 
across  the  heater.  A  fine  adjustment  was  provided  "by  a  one  ohm  variable 
resistor  in  series  with  the  heater.   The  variacs  and  the  variable  resis- 
tor were  mounted  on  a  panel  together  with  an  ammeter  and  voltmeter,  and 
appropriate  switches,  fuses  and  indicator  lights. 

The  power  supplied  to  the  test  cylinder  heater  was  measured  con- 
tinuously by  an  external  Weston  Model  k^2   double -current  range  wattmeter. 
The  meter  was  connected  in  such  a  manner  that  the  meter  indication 
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included  only  its  potential  circuit  losses.  On  low  scale,  used  to  measure 
power  below  30  watts,  the  internal  resistance  of  the  potential  circuit  was 
4083  ohms,  as  measured  by  a  high  precision  pdrtable  wheatstone  bridge.  On 
high  scale,  used  to  measure  power  above  30  watts,  this  resistance  was  8l80 
ohms.   The  resistance  of  the  heater  was  90. 07  ohms;  the  resistance  of  the 
heater  leads  was  0.17  ohms.   Thus,  with  the  meter  on  low  scale,  the  po- 
tential and  lead  circuit  loss  was  2.l6  per  cent  of  the  wattmeter  indication; 
on  high  scale  the  loss  was  1.09  per  cent  of  the  wattmeter  indication. 

2-d.   Temperature  measuring  apparatus 

The  eight  thermocouples  mounted  in  the  test  cylinder  walls  and  the 
two  additional  thermocouples  used  to  measure  ambient  temperature  in  the 
test  room  were  connected  through  a  Minneapolis -Honeywell  thermocouple 
selector  switch  to  a  common  cold  junction  and  a  portable  Rubicon  precision 
potentiometer.  The  cold  junction  was  placed  in  a  thermos  bottle  filled 
with  water  and  ice. 

The  differential  thermocouples  of  the  two  heat  meters  at  the  test 
cylinder  ends  were  connected  to  the  Rubicon  precision  potentiometer 
through  a  double -pole,  double -throw  knife  switch.   The  thermocouple 
selector  switch  was  put  on  open  circuit  when  it  was  desired  to  obtain 
the  heat  meter  temperature  differential  at  each  end  of  the  test  cylinder. 

Standard  tables  were  used  to  convert  the  potentiometer  readings  to 
temperatures.   It  was  considered  that  errors  associated  with  determining 
absolute  temperatures  and  temperature  differences  with  this  system  were 
within  0.1°F. 
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CHAPTER  IV 
EXPERIMENTAL  PROCEDURE 

The  experimental  program  was  carried  out  in  two  phases.   The  first 
phase,  called  the  calibration  phase,  was  a  series  of  tests  performed  with 
no  vibration  in  order  to  determine  the  heat  loss  by  conduction,  Q^,  to 
the  test  cylinder  supports.  The  second  phase,  called  the  vibration  phase, 
involved  the  collection  of  heat  transfer  data  with  the  test  cylinder  in 
horizontal  vibration. 

The  calibration  phase  was  carried  out  with  the  test  stand  and  test 
cylinder  in  the  horizontal  position  as  shown  in  the  photograph,  Fig.  10. 
The  test  section  heat  input  was  set  at  integral  multiples  of  5  watts 
from  5  to  50  watts .   Two  additional  tests  were  conducted  at  wattmeter 
settings  of  60  and  70  watts .   The  rated  accuracy  of  the  wattmeter  was 
0.5  per  cent  of  full  scale  deflection;  hence,  the  maximum  possible  error 
of  the  wattmeter  reading  occurred  at  the  5  watt  level  and  was  3' 5  Per 
cent.   The  possible  error  at  the  70  watt  level  was  0.5^  per  cent. 

Since  the  heat  input  to  the  test  cylinder  was  kept  constant  during 
a  given  test,  the  steady  state  temperatures  of  the  test  cylinder  were 
approached  exponentially  in  time.   The  criterion  established  for  terminat- 
ing the  calibration  phase  tests  was  a  maximum  change  in  any  individual 
test  cylinder  thermocouple  reading  of  0.2°F  and  a  maximum  change  in  heat 
meter  temperature  differential,  At^™,  of  0.1  F  over  a  period  of  one- 
half  hour.   In  practice,  a  steady  state  temperature  was  generally  achieved 
with  less  variation  than  the  above  requirements.   The  average  calibration 
test  run  extended  for  a  period  of  approximately  four  hours . 
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The   ambient  temperature,    t   ,   was   given  by  a  thermocouple  placed 
about  15"  below  the  test  section.      A  second  thermocouple,    suspended  in 
the  test  room  about   3  ft  above  and  behind  the  test  cylinder  was  monitored, 
and  generally  read  within  0.2°F  of  the   thermocouple  below  the  test  sec- 
tion. 

The   surface   temperature  of  the  test  section  was   defined  as   the 
average  of  the  six  temperature   readings   from  the   six  thermocouples   in  the 
test  section  locations   indicated  in  Chapter  III,    Section  2-b  .     As   expect- 
ed,   the   axial  distribution  of  the   individual  test  section  temperature 
readings   showed  a  variation  dependent  on  heater  power  level.      The  variation 
between  the  center  and  end  individual  test  section  temperatures   ranged 
from  k.O  F  at  the   5  watt  level  to  ^-0.3°F  at  the   70  watt  level.      The  tem- 
perature potential,    At,   based  on  average  surface   temperature,   t   ,    from 
the   six  test  section  thermocouple   readings  was   25.6  F  for  the   5  watt  level 
and  273.9  F  for  the   70  watt  level.      Thus,    considering  the  worst  condition   - 
that  the   average  tg  was   in  error  by  one-half  of  the   axial  temperature 
range    (2°F  for  the   5  watt  case   and  20.15  F  for  the   70  watt  case)   -  the 
percentage  error  in  At  would  be   J.Ql  per  cent  at  the    5  watt  level  and 
7.38  per  cent  at  the   70  watt  level.      These  errors,   of  course,    represent 
the  unlikely  extreme;    in  actual  test,    the  measured  average  test  cylinder 
surface  temperature  was   estimated  to  be  within  2  per  cent  of  the  true 
mean  surface   temperature . 

The  temperature  differential  across  the  heat  meters  at  the  test 
cylinder  supports  was  provided  by  the  differential  thermocouples  (described 
in  Chapter  III,  Section  2-*b).   The  temperature  level  from  which  the  tem- 
perature differential  was  measured  was  provided  by  the  thermocouple 
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immediately  beneath  each  heat  meter  in  the  test  cylinder  surface    (at  Sec- 
tions A-A  and  G-G,   Fig.    5)°      The   left  and  right  heat  meter  temperature   dif- 
ferentials were   always  within  2  F  of  each  other.      This  variation  was   due 
to   a  small  difference   in  heat  flow  across   the  left  and  right  test  cylinder 
heat  meters.      In  practice,    the   average  of  the   two  heat  meter  temperature 
differentials  was   taken  as   the  test  cylinder  end  temperature  difference, 

At        .As   stated  in  Chapter  III,    Section  2-d,    it  was   considered  that 
hm 

the  errors   associated  with  determining  absolute   temperatures   and  tempera- 
ture differences  with  this   system  were  within  ±  0.1  F. 

The   results  of  twelve   calibration  phase   test  runs   are   given  in 
Table   1  and  plotted  in  Fig.    13.      The  method  used  in  calculating  these 
results   is  presented  in  Chapter  VI. 

The  vibration  phase  of  the  experimental  program  was  oriented  toward 
the  major  objective  of  the  entire   investigation;    namely,    the   investigation 
of  the   influence  of  horizontal  mechanically  induced  simple  harmonic  vibra- 
tions upon  the  heat-transfer  rate  by  free   convection  from  a  heated  hori- 
zontal cylinder  to   air.      Specifically,    the  tests  were  so   arranged  as   to 
provide  data  which  could  be   compared  directly  with  the  sound  method  data 
obtained  by  Fand  and  Kaye   as  explained  in  Chapter  I. 

The   three  principal  experimental  variables  were   as   follows: 

a.  temperature  potential,  ^t 

b.  intensity  of  vibration,  af 

c.  heat  input  to  the  heat  transfer  cylinder  test  section,  Q 

In  order  to  provide  a  direct  comparison  with  the  sound  method  experiments 
of  Fand  and  Kaye  in  a  small  number  of  test  runs,  it  was  decided  to  restrict 
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the  independent  variable  At  to  one  value,  namely,  At  =  150  F.  Since, 
in  practice,  it  would  "be  extremely  difficult  to  obtain  an  exact  steady 
state  value  of  /\.t  =   150  F  for  each  test,  the  allowable  variation  from 
this  temperature  potential  was  established  as  ±  1.0  F.   It  was  also  nec- 
essary to  vary  af  within  a  range  of  about  0.2  to  1.1  ft/sec  to  cover  the 
same  range  of  vibration  intensity  as  was  done  in  the  sound  method  experi- 
ments . 

The  vibration  tests  we're  all  made  at  the  same  frequency,  f  =  104-  cps, 
the  resonant  frequency  of  the  mechanical  system.  At  the  start  of  each 
test,  the  vibration  exciter  was  tuned  to  this  resonant  frequency.   Tuning 
was  achieved  by  applying  the  principle  that  for  a  coupled  electro- 
mechanical system  with  a  constant  electrical  energy  input,  the  current 
becomes  a  minimum  at  resonance .  An  additional  aid  in  determining  reson- 
ance was  the  fact  that  the  maximum  apparatus  noise  level,  while  not  large, 
coincided  with  the  resonant  condition. 

After  obtaining  vibration  at  resonance,  the  desired  value  of  af  was 
acquired  by  adjusting  the  output  level  of  either  the  audio -oscillator  or 
its  amplifier  or  both.  Also,  at  lower  vibration  intensities  it  was  found 
convenient  to  limit  the  value  of  the  direct  current  output  to  the  exciter 
field  coil  to  an  arbitrary  value  less  than  the  1.75  ampere  rating.   This 
procedure,  in  effect,  added  sensitivity  to  the  gain  controls  of  the 
audio -oscillator  and  amplifier  and  thereby  aided  in  establishing  the  de- 
sired amplitude.   These  gain  controls  were  adjusted  in  trial  and  error 
fashion  until  the  double  amplitude  of  vibration  was  at  the  value  required 
to  give  the  proper  af  value  for  the  test  run.   This  procedure,  initially 
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tedious,  was  usually  accomplished  within  10  minutes  after  some  experience 
had  been  acquired.   The  vibration  amplitude  was  checked  at  5  minute  inter- 
vals; minor  correction  was  occasionally  necessary. 

The  requirement  of  producing  a  steady  state  value  of  A  t  within  1  F 
of  150  F  was  more  difficult  to  achieve.  Because  the  steady  state  heat 
input  to  the  test  cylinder  with  vibration  was  unknown,  its  value  had  to 
be  determined  in  each  test  run.  The  power  supplied  to  the  test  section 
which  resulted  in  a  150  F  temperature  potential  in  the  absence  of  vibra- 
tions was  determined  to  be  about  3^-  watts  from  the  results  of  the  calibra- 
tion tests.  It  was  assumed  that  the  power  level  needed  to  produce  the 
same  temperature  potential  with  the  test  cylinder  in  vibration  was  greater, 
The  test  cylinder  was  heated  for  about  two  hours  prior  to  vibration  at  a 
power  level  somewhat  greater  than  3^  watts .  After  the  desired  vibration 
intensity  was  established,  a  complete  set  of  temperature  readings  was 
taken  at  10  minute  intervals .   The  temperature  potential  existing  at  the 
time  of  temperature  measurement  was  then  quickly  computed,  and  the  power 
level  was  adjusted  to  bring  At  closer  to  150  F.   This  process  was  re- 
peated until  no  change  in  A  t  greater  than  0ol°F  in  a  10  minute  period 
was  detected.  At  this  time  the  watt  meter  reading  was  recorded  as  the 
steady  state  value,  and  the  cathetometer  readings  of  vibration  amplitude 
were  recorded. 

The  duration  of  the  individual  vibration  tests  varied  between  two 
and  three  hours .   The  results  of  the  nine  tests  with  vibration  are  re- 
corded in  Table  2  and  plotted  in  Fig.  1^-.   The  method  used  in  calculating 
these  results  is  presented  in  Chapter  V. 
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CHAPTER  V 
CALCULATION  PROCEDURE  AND  RESULTS 

This  chapter  contains  a  description  of  the  methods  and  basic  calcula- 
tion procedure  used  to  evaluate  and  interpret  the  experimental  data.,  and 
a  description  of  the  results  obtained. 

The  calibration  phase  tests  of  the  experimental  procedure  provided 
data  with  which  the  axial  heat  losses  from  the  test  cylinder  could  be  de- 
termined as  a  function  of  average  heat  meter  temperature  differential,,  At,  . 
During  the  tests  in  the  absence  of  vibration,,  the  heat  input  to  the  test 
section,  QL.,   was  dissipated  to  the  system  environment  by  the  following  modes 
of  heat  transmission: 

a.  Heat  loss  from  the  test  section  by  radiation,  Q~ 

b.  Heat  loss  from  the  test  section  by  convection  (in  the  absence 

of  vibration),  Q 

"   CO 

c.  Heat  loss  to  the  test  cylinder  supports,  Qj 

d.  Heat  loss  by  convection  from  the  end  terminal  plates  of  the 
test  cylinder  and  heat  loss  by  conduction  through  the  test 
cylinder  lead  wires  (in  the  absence  of  vibration),  Q^q?    called 
the  unaccountable  heat  loss. 

The  test  section  heat  losses  by  radiation  and  convection,  Qr  and 
Qco,  were  subject  to  accurate  calculation  as  will  be  described  below. 
Calculation  of  the  heat  loss  to  the  test  cylinder  supports,  Qd,  sould  not 
be  accurately  calculated,  however,  because:   (l)  the  thermal  conductivity 
of  the  heat  meter  Teflon  ring  was  not  known  within  the  required  accuracy,, 
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and   (2)   the   contact  coefficient  of  the   copper-Teflon  heat  meter  junctions    — 
a  function  of  test  cylinder  support  clamping  pressure,    surface   roughness 
and  the  presence  of  oxide   films    —  could  not  be  determined.      Further,    cal- 
culations  for  Q-, ,   using  the  "best  available  value  of  thermal  conductivity 
for  Teflon,    gave   results   that  did  not  agree  well  with  the  empirical  deter- 
initiation  of  Q,-,   as   described  below.      Similarly,    Q^0  could  not  be   calculated 
because  of  the   complicated  geometry  of  the  outside   surface  of  the  Teflon 
plates . 

The   sum, (0^  +   Quo)*    could  be  determined,    however,    through  application 
of  the   following  heat  flow  continuity  equation: 

Hi+    %x>      =     %  -   (Or  +    Qco)  (6) 

In  the  calibration  phase  tests,  (Q^  +  Qnn)  "was  determined  as  a  function 
of  At^.   The  results  of  these  tests  are  plotted  in  Fig.  13.   The  origi- 
nal experimental  data  for  the  calibration  tests  and  the  associated  cal- 
culated quantities  are  listed  in  Table  1.   The  experimentally  determined 
curve  of  (0^  +  Quo)  versus  .A^hm  ^s  very  nearly  a  straight  line  which 
indicates  that  the  predominate  heat  loss  at  the  test  cylinder  ends  was  by 
conduction,  i.e.,  Q^  +  Quo"   This  was-  the  expected  result  since  the  out- 
side surface  temperature  of  the  Teflon  terminal  plates  was  much  less  than 
the  temperature  of  the  test  cylinder.  Also,  the  heat  loss  by  conduction 
through  the  lead  wires  was  negligible  compared  to  Q, . 

The  heat  transferred  by  convection  in  the  presence  of  vibrations,, 
Qcv,  was  determined  from  the  heat  continuity  equation: 

Qcv  =  %  -   Qr  -  (Qd+  Quv)  (7) 
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where  Quv  is  the  unaccountable  heat  loss  in  the  presence  of  vibration.  It 
is  recognized  that  Quv  was  somewhat  greater  than  Quo  due  to  the  increase 
in  convective  heat  transfer  from  the  terminal  plates  in  the  presence  of 
vibrations;  however,  this  difference,  (Quv  -  Quo)  was  assumed  negligible 
since  the  unaccountable  heat  loss  itself  was  small.  Accordingly,  in  the 
tests  with  vibration,  Qcy  was  computed  by: 

Qcv  =  ^  -  ^  -  (Qd  +  Qu0)  (8) 

where    (Qq  +   Q^q)  was   determined  from  Fig.    13 . 

The  heat  transferred  by  convection  from  the  test  section  in  the  ab- 
sence of  vibration,  Qco,  was  determined  in  the  following  manner.  The  de- 
fining equation  of  the  heat-transfer  coefficient  in  the  absence  of  vibra- 
tion,   hQ,    is: 


Q 


co  o 


h_      =     3.1H2 — ,  Btu/hr-ft   -°F  (9'0 

°  A  At 


where 


2 
A  -     test  section  heat  transfer  area,    0.0993  ft 

Q  =     heat  transferred  by  convection,   watts 

At        =     temperature   potential,      F 

3.412  is   the   conversion  factor  from  watts   to  Btu/hr 

(9) 
The  experimentally  determined  equation         for  hQ  of  a  cylinder  of 


diameter  D     is 


At     XA 

h0      =     0.225   (~d~")  ,    Btu/hr-ft2-°F  (10) 

o 
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Equation   (10)  was   the   result  of  a  heat-transfer  study  with  a  cylinder  of 
identical  diameter   (7/8"),    length  and  end  support  method  and  surface  thermal 
emissivity  as   the  test  cylinder  used  in  this   project. 

Equating  Equations    (9)   and  (10)  with  DQ     =     0.0729  ft  and  solving  for 
Qco  gives:  ^ 

Qco  =  0.01^3  A   t     ,  watts  (11) 

The  maximum  experimental  error  in  the  determination  of  hQ  by  Equa- 
tion (l0)  was  reported  to  "be  2.6  per  cent.     The  maximum  possible  experi- 
mental error  in  the  calculated  value  of  Qc0  was  determined  by  logarithmic 
differentiation  of  Equation  (9): 

dh       dQ       dA     d(^t) 

O  CO  .    . 

h0         %0  -    a  ■  "XT  ll2) 

Equation  (12)  was  written  in  finite  difference  form;  hence ,  the  maximum 
possible  percentage  error  in  Q   was  the  sum  of  the  maximum  possible  per- 
centage errors  of  the  quantities  involved  in  its  calculation.   The  possible 
error  in  h^   was  2.6  per  cent;  the  possible  error  in  the  calculation  of 
surface  area,  A,  was  less  than  0.2  per  cent.   The  maximum  error  in  ^t  was 
estimated  to  be  2  per  cent  (see  Chapter  IV).  Thus,  the  maximum  possible 
error  in  the  value  of  Q,c0  was  k.8   per  cent. 

The  heat  transferred  by  radiation  from  the  test  section,  Qr,  was 
independent  of  the  presence  of  vibrations.   Q  ,  as  a  function  of  At,  was 
determined  from  the  graph,  Fig.  12,  based  on  the  equation: 

Qr     =     0.2931  CTe  A 


(tb  +        t)      -  tb        ,   watts        (13) 
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where 

C~  =     Stefan -Bo It zman  constant,  0.1713xlO~8  Btu/ft2-hr-°R^ 

e       =  surface  thermal  emissivity  of  the  test  section,  0.0^5 

p 
A  =     test  section  heat  transfer  area,    0.0993  ft 

t  =     mean  ambient  temperature,    5^0  R 

^  t     =  temperature  potential,  F 

0.2931     -     conversion  factor  from  Btu/hr  to  watts 

The  errors   associated  with  determining  Q  ,   which  were   introduced  by  the 
fact  that  the   actual  ambient  temperature  was   in  most  cases   slightly  dif- 
ferent from  the  mean  ambient  temperature,   were   less   than  k-  per  cent. 
These  errors  were  negligible  because  Q     itself  was   small  compared  to  the 
total  heat  transferred  from  the  test  section. 

The  heat  transfer  coefficient  in  the  presence  of  vibration,    h   ,   was 
calculated  for  each  test  of  the  vibration  phase  by  using  the   heat  trans- 
ferred by  convection,   Qcv>    for  the   test  as   determined  by  Equation   (8). 

Thus: 

Oct 
\    =  3.412  Y&t     >   watts  (li^ 

The  maximum  possible  experimental  error  in  the  calculated  value  of 
h  was  determined  by  writing  the  logarithmic  differential  equation  of  (ik) 
in  finite  difference  form  as  in  the  case  for  Qco  above.   The  maximum  pos- 
sible error  involved  in  obtaining  Qcv  was  3-0  per  cent,  as  obtained  by 
adding  the  maximum  errors  possible  in  Q^  and  (Qd  +  Quo).  Errors  in  Qr 
were  less  than  0.1  per  cent  and  were  neglected.  Again  the  errors  in  sur- 
face area,  A,  and   ^  t  were  taken  as  0.2  and  2  per  cent,  respectively. 
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Thus,    the  maximum  possible  experimental  error  in  the  determined  value  of 
h     was    5.2  per  cent. 

The  experimental  data  and  the  calculated  results  of  nine   tests  with 
vibration  are   recorded  in  Table  2.      The   range  of  vibration  intensity  em- 
ployed was   from  0.213  to   1.109  ft/sec;    the  temperature  potential  in  all 
these  tests  was  within  ±  1  F  of  150°F.      The  value  of  hy  varied  from 
1.803  Btu/hr-ft2-°F   at  the   lowest   af  to   4-323  Btu/hr-ft2-°F  at  the  highest 
af .      The   value   of  hQ   at    ^  t      =     150°F  was    computed  to  be   1.717  Btu/ 
hr-ft   -°F  by  Equation   (l0).      Thus,    percentage   increases   in  the  heat-transfer 
coefficient  from  5.0  to   151.8  per  cent  were   achieved. 

In  order  to  provide   a  direct  comparison  of  these   results  with  the 
sound  method  investigation  of  Fand  and  Kaye        '    (discussed  in  Chapter  I), 
the  variation  of  h^  with  af  at    ^  t     =     150  F  for  a  stationary  horizontal 
cylinder  in  the  presence  of  a  high  intensity  sound  field  of  frequency 
f     =     1496  cps  was  obtained  from  Fig.    6l  of  Reference   10  by  cross -plotting. 
The  data  points  of  the  vibration  tests  of  the  present   investigation  were 
then  plotted  on  the   same   graph,    Fig.    Ik.      The  maximum  deviation  of  any 
single  data  point  from  the   sound  method  curve   in  Fig.    lk  is   k.6  per  cent 
as   compared  to   the  maximum  possible  experimental  error  of   5.2  per  cent  in 
obtaining  the  data  points .     Fand  and  Kaye   arrived  at  a  correlation  equa- 
tion which  expressed  the  dependence  of  h     on    /\  t  and  af  with  a  high  degree 
of  consistency  above   an  af  of  about  0.7  for  sound  frequencies   less   than 
2000  cps.      For  convenience,    this   correlation  equation  is   restated  here: 

1/3 


h       =     0.722 
v 


Z\t    (af)2 


(1) 
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The  maximum  deviation  of  h  obtained  in  this  investigation  for  the  three 
values  of  af  greater  than  0.7  from  the  computed  value  h  given  by  Equa- 
tion (l)  was  7*3  VeT   cent.  The  average  deviation  was  2.1  per  cent. 

A  comparison  of  the  results  of  this  investigation  with  the  results 

of  the  study  of  heat-transfer  rate  variation  with  vertical  mechanical 

(8) 
vibrations  (discussed  in  Chapter  I)  was   also  made.      A  curve  of  h     ver- 

sus  af  for    A  "t     =     150  F  was  obtained  from  Fig.    12  of  Reference   8  and  is 
plotted  in  Fig.    15  of  the  present  report.      The  data  points  of  the  vibra- 
tion phase  of  this   investigation  were  plotted  on  Fig.    15  to   facilitate   com- 
parison.     Inspection  of  Fig.    15  revealed  that  the   trend  of  h^  increase 
relative   to   af  was   quite   similar  for  the  two   investigations.      In  particu- 
lar,   at  values  of  af  above   about  0.3  ft/sec,    the  value  of  h     is   somewhat 
less   for  the   horizontal  vibration  case,   while  the  slopes   correspond  closely. 
This   result  will  be  discussed  further  in  Chapter  VI. 

It  was   considered  desirable  to  obtain  some  visual  evidence  of  the 
type  of  boundary-layer  flow  that  existed  around  the  heated  test  cylinder 
during  vibration.     Accordingly,    a  crude   flow-visualization  examination  was 
performed,   using  a  lighted  cigarette   held  beneath  the   cylinder  as   a  smoke 
source .      The   results  of  this  visual  observation  clearly  indicated  that 
laminar  flow  persisted  in  the   region  near  the   lower  stagnation  point 
throughout  a  range  of  af  from  0  to   about  1.1  ft/sec.     Also,    some  evidence 
of  vortex  formation  above  the   cylinder  was  noted;    this  vortex  flow  is 
characteristic  of  thermoacoustic  streaming  as   observed  by  Fand  and  Kaye 
in  their  sound  studies.      Although  the   results  of  the  elementary  visualiza- 
tion technique  employed  here   require   further  substantiation,    the  study 
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did  indicate  conclusively  that  the  fluid  flow  mechanism  in  the  region  of 
the  lower  stagnation  point  of  the  cylinder  was  not  vib rationally  induced 
turbulence  as  reported  for  the  vertical  vibrations  case   ' 
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CHAPTER  VI 
DISCUSSION  OF  RESULTS 

The  results  of  the  present  experimental  investigation  reaffirm  the 
conclusions  of  previous   investigations   that  the   influence  of  vibrations 
on  heat  transfer  is   significant.      The  plots  of  h     versus   af   (Figs.    1^  and 
15)   clearly  indicate   the   general  effect  of  the   intensity  of  vibration  on 
the  heat-transfer  rate   from  a  heated  horizontal  cylinder  in  horizontal 
transverse  vibration  in  air.      It  was  noted  that  the   change   in  h.  at  vibra- 
tion intensities   less   than  af  -  0.3  ft/sec   is  not   appreciable.     Above  this 
value  of  vibration  intensity,    the  heat-transfer  coefficient  increases 
rapidly  with  increasing  af .      The  value  of  af  =  0.3  ft/sec,   where  vibra- 
tions  take  on  major  significance  with  regard  to  the  heat-transfer  rate 
is   called  the   critical  intensity.      The   critical  intensity,    af  =  0.3  ft /sec, 
observed  in  this   investigation  corresponds  very  nearly  to   the   critical 
intensities   reported  for  the  sound  method   (Fig,    Ik)   and  also   in  the  ver- 
tical mechanical  vibrations  method   (Fig.    15)» 

The   close   agreement  of  the   quantitative  heat  transfer  results   and 
the   flow  visualization  observations  of  this    investigation  with  those  of 
the   sound  method  for  comparable  vibration  vector  orientation  provides 
evidence  which  indicates   that  the  physical  mechanism  of  the  heat  transfer 
for  the  two   cases   is   identical.      However,    the   quantitative  results,   when 
compared  to   the  vertical  vibrations  method   (Fig.    15) j    show  that  the 
slopes  of  h     versus   af  are  very  nearly  the  same   for  these  two  cases.      The 
fact  that  the   trend  of  the  data  obtained  in  this   investigation  is   similar 
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to  the  results  of  the  vertical  vibrations  method  (Fig.  15)  is  not  con- 
sidered evidence  that  the  two  modes  of  vibration  are  equivalent  with 
respect  to  heat  transfer;  for,  as  has  been  explained  in  Chapter  V,  the 
mechanism  of  interaction  between  vibrations  and  free  convection  in  the 
case  of  vertical  vibration  is  vibrationally  induced  turbulence,  whereas 
there  is  evidence  that  the  mechanism  of  interaction  for  both  the  sound 
method  and  the  present  investigation  is  thermoacoustic  streaming. 

The  slope  of  h  versus  af  for  af  greater  than  about  0.9  for  the 
v 

sound  method  is  2/3  as  Equation  (l)  indicates.   The  corresponding  slope 

(9) 
for  the  vertical  vibrations  method  has  been  reported  to  be  1   .  Addi- 
tional heat  transfer  data  is  required  for  values  of  af  greater  than  0.9 
to  determine  conclusively  the  value  of  this  slope  for  the  horizontal 
vibrations  method. 

The  ratio  of  amplitude  of  vibration  to  cylinder  diameter,  a/D , 
is  a  controlling  variable  in  all  fluid-mechanical  phenomena  involving 
vibrations .   The  a/D  ratio  for  this  investigation  at  a  typical  value 
of  af  =  0.9  ft/sec  and  f  =  10U  cps  was  0.119;  the  a/DQ  ratio  for  the 
sound  method  of  Fand  and  Kaye  with  af  =  0.9  ft/sec  at  f  =  1^96  cps  was 
9.63  x  10"3.   Thus,  the  a/DQ  ratios  differed  by  a  factor  of  12. 3.   In 
view  of  the  close  quantitative  agreement  of  the  results  of  the  two 
methods  of  investigation  (Fig.  lh) ,    it  appears  that  the  influence  of 
a/DQ  does  not  have  a  predominant  effect,  at  least  for  the  range  of  a/DQ 
ratios  indicated  here . 
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CHAPTER  VII 
CONCLUSIONS 

The  experimental  portion  of  this   investigation  of  the   influence  of 
horizontal  mechanically  induced  vibrations  on  the  heat-transfer  rate   from 
a  horizontal  cylinder  to  air  consisted  of  a  total  of  21  heat  transfer 
test  runs.      Twelve  of  these  test  runs  provided  calibration  data  on  axial 
heat  loss   from  the   heat  transfer  test  cylinder;    nine   test  runs  were   con- 
ducted with  the   cylinder  in  horizontal  vibration.      The   ranges   of  the 
principal  experimental  variables  were: 

temperature  potential,   At:  150°F  +.   1  F 

frequency  of  vibration,    f :  10U  cps 

amplitude  of  vibration,    a:  0.025  to   0.128" 

intensity  of  vibration,    af:  0.21  to   1.11  ft /sec 

The   results   indicate   that  the   influence  of  horizontal  vibrations 
upon  the  heat-transfer  rate  by  free  convection  to  air  becomes   appreciable 
when  the   intensity  of  vibration,   defined  as   the  product  af,   exceeds   the 
critical  value  of  0.3  ft /sec .      This   critical  value   is   approximately  equal 
to   the   critical  sound  vibration  intensity  determined  in  Reference   10  and 
the  critical  value  of  af  for  vertical  mechanical  vibrations   determined  in 
Reference   9-      The   flow  near  the  horizontal  cylinder  in  the  presence  of 
intense  mechanically  induced  transverse  horizontal  vibrations  was  ob- 
served to  be   laminar  near  the   lower  stagnation  point  using  smoke   as   the 
indicating  medium;    this   type   of  flow  is  one   characteristic  of  the   flow 
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near  the  lower  stagnation  point  of  a  cylinder  in  the  presence  of  intense 
acoustically  induced  transverse  horizontal  vibrations.   The  observed 
laminar  flow  near  the  lower  stagnation  point  was  quite  different  from  the 
vibrationally  induced  turbulent  flow  characteristic  of  horizontal  cylin- 
ders in  transverse  vertical  vibration. 

Based  on  the  quantitative  heat  transfer  results  and  the  elementary 
flow  study  of  this  investigation,  it  appears  that  the  physical  mechanism 
of  heat  transfer  for  the  sound  and  mechanical  vibrations  methods,  both 
with  geometry  as  indicated  in  Fig.  1-a,  is  identical.  Also,  the  com- 
parison of  these  two  methods  indicates  that  the  amplitude  of  vibration 
to  cylinder  diameter  ratio,  a/DQ,  does  not  have  a  predominant  effect  for 
the  range  of  a/D  ratios  considered. 
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CHAPTER  VIII 
RECOMMENDATIONS 

This  chapter  presents  recommendations  for  future  investigations  in 
the  field  of  heat  transfer  from  vibrating  cylinders  based  on  the  results 
and  experience  gained  from  the  present  project.   Six  specific  suggestions 
are  made;  the  first  five  of  these  relate  to  experimental  work  that  can  be 
conducted  with  the  apparatus  designed  and  constructed  for  this  project. 
Indeed,  this  apparatus  was  designed  to  include  features  whereby  adaptation 
to  other  investigations  will  require  minor  or  no  modification.   The  list 
of  suggested  investigations  with  a  few  relevant  remarks  is  given  below. 

(1)  Extension  of  the  range  of  vibration  intensity  to  values  of  af 
greater  than  1.1  ft/sec.  This  is  necessary,  as  indicated  in  Chapter  VI, 
in  order  to  determine  accurately  the  slope  of  h  versus  af  in  the  higher 
af  range.  The  existing  apparatus  has  been  tested  at  af  =  1.3  ft/sec 
with  about  35  lbs  exciter  driving  force.  It  is  considered  that  an  af  of 
at  least  1.5  ft/sec  can  be  achieved  at  the  maximum  exciter  output  of  50 
lbs  force , 

(2)  Flow  visualization  studies.   It  would  be  of  particular  value  to 
obtain  photographic  evidence  of  the  flow  pattern  around  the  test  cylinder. 

(3)  Change  in  vibration  direction  relative  to  the  direction  of  the 
force  of  gravity.   Since  the  flow  character  around  the  cylinder  has  been 
shown  to  be  different  for  the  vertical  and  horizontal  vibration  cases,  it 
would  be  of  value  to  investigate  vibration  angles  between  these  two  ex- 
tremes.  This  can  be  done  with  the  existing  apparatus  as  described  in 
Chapter  III,  Section  1. 
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(k)     An  investigation  of  velocity  and  temperature  distributions 
near  the  cylinder.   It  is  considered  that  hot-wire  anemometry  techniques 
could  be  used  to  determine  detailed  information  concerning  the  heat 
transfer  mechanism  near  the  cylinder  by  an  investigation  of  this  type. 

(5)  Different  cylinder  diameters.   The  quantitative  determination 
of  the  influence  of  cylinder  diameter  is  important.   New  test  cylinders 
and  some  modification  of  the  test  cylinder  yoke  would  be  required. 

(6)  Analysis.   It  is  considered  that,  on  the  basis  of  the  data 
provided  herein  and  from  other  investigations  in  the  field,  an  analysis 
should  be  attempted  which  would  provide  a  mathematical  foundation  for 
the  phenomena  observed. 
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APPENDIX 


TO. 

APPENDIX  A 
VIBRATION  ANALYSIS 

In  the   analysis   of  the   free-free  beam  and  coupled  center  load  reso- 
nant vih ration  syste,    it  was   assumed  that  the  exciter  moving  element, 
the   test  cylinder  and  its   supporting  yoke   and  the   hardware  used  in  coup- 
ling the  exciter  moving  element  to   the  beam  acted  as   a  concentrated  mass 
of  magnitude  m  at  the  beam  center.      The  beam  had  a  length  2.J[   and  a 
total  mass  mu .      The   coordinates  used  in  the   analysis   are   as   shown  below. 


x 


-X 


x  =   0 


x   =  2, 


/ 


l\ 

■t 

^^ 

"^•^                                      ^"^ 

V 

x  = 

=    X 

n 

->-  x 


The  general  equation  for  the  transverse  vibration  of  a  prisimatical 

(13) 
beam,  as  presented  by  Timoshenko     is: 


dx^ 


<a  a2: 


go  at2 


(15) 
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which  may  "be  "written: 


j  2  Ele        ,k 

d  y  &o     dH; 


+ 


dt*         aY     a4* 


=    o  (16) 


where 


E        =     modulus  of  elasticity,    lb/in 

k 
I        =     beam  cross-sectional  moment  of  inertia,    in 

o 
g   =  gravitational  constant,  386  in/sec 


q       =     beam  weight  density,    lb/in 


Considering  only  the   first  mode  of  vibration  as  was   justified  in 
Chapter  III,    Section  1-b,    the   solution  to  Equation   (l6)   has   the   form 
y     =     X  sincOt,  ,   where   sinoot       gives   the  periodic  variation  of  the 
beam  deflection  along  its   length  according  to   the   function  X  =  f(x). 
Substituting  this   solution   into  Equation   (l6)  yields: 

4"      <A     -     0  (17) 

dx 

where  o(   is  defined  by: 

0(.U    -    ^-  (18) 

ET8o 
Particular  solutions  of  Equation   (l8)   involve   the   functions   sin^y    , 

cos  0*  X      ,    sinhoCX      and  coshp'y      .      Therefore,   the  solution  to  Equation   (l6) 

can  be  conveniently  expressed  as: 


y     =     X  sin 


LdZ  =  C-j_   (cos<XX    +   cosh<AY    )  +   C2   (cospty.   -  cosh^y.   ) 

+   C~    (sinfltX  +  sinhpty    )  +   Ck    (sin</y    -  sinhyx   )   sinuoiT 

(19) 
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Application  of  the  end  boundary  conditions,    lee.,    the  bending 

d  y  cj3v  a 

moment     EI     — o     and  shear  force     EI  — *-     vanish  at  x  =  0  and  x  =  Zyf  , 
dx^  dx3 

reveals   that  Cp   =  Cr    =  0.      Similarly,    the  mid-span  boundary  condition 


dy 
of  zero   slope,     —     =     0  at  x 
dx 


-I  . 


gives 


cos 


dA    +  cosh^y 


sinof/    -  sinhtf/  3 


(20) 


An  additional  condition  at  beam  center  is  obtained  from  the  force - 
inertia  balance  at  x  =  a     : 


(a)  Exciter  force    (up): 

(b)  Shear  force    (up): 


(c)      Inertia  load  of  m     (down) 


Thus ,    at  x  =    / 


F     =     F0   sin 60 "t 


2EI 


dx3 


asy 

-  m  o 

c     dt^ 


-t  d  y 

F     sint^T    +   2EI  — * 

o  -,   ~> 

dx 


m 


d2y 
c     dt2 


(21) 


Substitution  of  Equations  (19)  and  (20)  into  Equation  (21)  results 
in  the  following: 


F  +   2EI  0(       Co 


-   cos 


\u>     c3 


ftjf  +   coshrt/  +    (sino</  +   sinh#/  ) 
sin«.4    +   sinho</   + 


<X>f  +   coshft<f 


sin«/  -   sinh 
(costf^    +   coshc*/    ) 


sintfyf    -   sinhrt^f 


(22) 
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Also,  from  the  definition  of  </      the  beam  mass,  m,  ,  may  be  conveniently 
introduced  into  Equation  (22)  as  follows: 


-*3  ■  ¥  ^  ■  w 


For  convenience,  the  first  and  second  terms  in  brackets  of  Equation  (22) 
are  called  Q  and  R.   Substituting  Equation  (23  into  (22)  gives: 

S   =  "^ (2^ 

Examination  of  Equations    (19)   and  (22)   shows   that  beam  displacement, 
y,    is   a  function  of  CU   and      MA  .     At  the  theoretical  resonance   condition, 
y  becomes   infinite.      However,    for  a  finite   frequency  and  a  beam  of  finite 
mass   and  dimensions,    Oil  cannot  approach  infinity.      Hence,   with  regard  to 
theory,    CU  must  become   infinite   at  resonance.      Therefore,    the  denominator 
of  the   right-hand  side  of  Equation   (24)  must  vanish  at  resonance,    and: 


"\ 


mc 


<*/    "  (25) 


Substitution  of  the  bracket  terms   Q  and  R  and  simplifying  produces   the 
following  frequency  equation  in  terms   of  beam  to   center  mass   ratio   at 
resonance: 


"b  OLA    (sechflU     +   costf/ 


m 
c 


(sintf/      +    cosotf   tanhrf/)  ^) 


Equation   (26)    is  plotted  for  various  mass   ratios  versus   values  of  Q(J(    in 
Fig.    6.      It  is  worthwhile   to  note   that  the   first  mode   solution  of  the 
simple   free-free  beam  with  no   added  center-span  mass   should  correspond  to 
the  value  of     (/J[        as  m^/m     approaches   infinity.      Reference  Ik  gives   the 
frequency  solution  for  this   case   as: 


EIgQ 

U)  =      22.4      -\j       2_  (27) 

AV(2i)4 

Substitution  for    p(     into  Equation   (27)   and  solving  for    o(Ji      gives   a 
value  of  0(a      -  2-37  at  resonance.      Fig.    6  clearly  indicates   that 
approaches   2-37  as   n^/n^   approaches   infinity. 

In  order  to  have  minimum  force   transmission  from  the  beam  to   the 
beam  supports,    these   supports   must  be   located  at  the  vibration  npdal 
points.      Solving  for  beam  deflection,    y,      by  substituting  Equation   (20) 
into    (19): 

costf^    +   coshpM      ,  ,     .     .  ,  ,~ 

sinetf    -   sinhd>      (C°S*X   +    cosh«*    )+(sln0()(    +    SlnhC?(^ 

C     sinfcOt    =     0  (28) 

At  the  nodes,    say  at  x  =  xn,   y  =  0.     Equating  Equation   (28)   to   zero   gives 


sin 


OlJfr        +   sinh  <*/_!" 

cos  (/J  +   cosh^-f 


cos  lA 


.1  ' 
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The  ratio  of  nodal  distance  from  beam  end  to  "beam  half  length,  x  / J    ,  is 

plotted  versus  0(J(      in  Fig.  J. 

The  maximum  "bending  moment  in  the  beam  will  occur  at  x  =  Jf     when 

the  center  amplitude,  a,  is  at  its  maximum  value.   This  will  take  place 

at  times  when  slnCOt   is  unity.   From  elementary  beam  theory,  the  bending 

d2 

moment   is   given  by  M  =  EI  — 75  ;    thus,    differentiating  Equation   (28)   twice: 

dx 


2 

M  =  EI  — J 

max  dx 


2  sinOM    sinh#>, 

-     ^X      c,   (sin^  -sinhtf/  )      (30) 


==/ 


The  maximum  stress,  ^""max  >    i-s  given  by  M   /i  times  h/2  where  h  is  the 
beam  height  in  the  direction  of  deflection.   Combining  this  with  Equations 
(2k   and  (28)  gives: 

Oi-'*£^j*-&       ^b    u2    +    U3  (31) 

mc 


where 


2 
I 

3 
u 


U       =     (^)3  sin  (Xyftanh^/  (32) 

U2     =     sintf/  +   cos  &J  t&nhtff  (33) 

Uo     =      tf/(sechflf/  +   cos  a./    )  (3*0 


-,,   U     and  U_   are  plotted  versus    Q(J      in  Fig.    8. 

The  maximum  amplitude   at  x  =/       is   given  by  Equation   (28)   for  times 


when  sin#)£    is  unity.     Upon  substitution  for  C    • 
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(35) 


F        1 

U 
3 

mc   CO2 

m       2 
c 

U3 

Equations    (31)   and  (35)   are  not  useful  at  resonance  because   the 

m, 

denominator  term  —  \J     +   U  vanishes  at  the  resonant  condition.   This 

m    ^    3 
c 
is  a  result  of  the  earlier  assumption  that  maximum  amplitude  goes  to 

infinity  at  resonance,  i.e.,  Co  approaches  infinity  at  resonance.   However, 

the  ratio  of  bending  moment  to  amplitude  may  be  obtained  by  dividing 

Equation  (31)  by  (35)>  and: 

O  max     E   h     1  ,  ., 

The  maximum  bending  stress  versus  amplitude  and  vibration  intensity,  af , 
for  steel  and  aluminum  beams  at  various  fXJf    values  and  at  f  =  100  cps  is 
plotted  in  Fig.  9« 

The  interpretation  of  this  analysis  and  final  beam  selection  is 
presented  in  Chapter  III,  Section  1-b .   The  curves  developed  by  means  of 
this  analysis  (Figs.  6,  7>  8  and  9)  "were  used  for  final  beam  design  ac- 
cording to  the  calculation  method  indicated  in  Appendix  B.   It  is  con- 
sidered that  this  analysis  served  as  a  remarkably  useful  tool  in  this 
project.   The  resonant  frequency  obtained  of  the  actual  beam  was  within 
1  l/2  per  cent  of  the  design  value. 
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APPENDIX  B 
RESONANT  BEAM  DESIGN  CALCULATIONS 

The   resonant   free-free  beam  design  calculations  based  on  the 
analysis  presented  in  Appendix  A  are   given  below.     Beam  w±c|bh,  b,   was 
initially  set  at  1.50"    so  that  existing  hardware   for  test  cylinder  sup- 
port could  be  employed.      The  material  selected,    "HY-100"   steel,    had  a 
weight  density  of  O.283   lb/in-^.      As   explained  in  Chapter  II,    Section  1-b, 
a  value  of  Q(J[   =2.25  was   selected  for  the   final  design  on  the  basis   of 
stress   and  beam  size   considerations.      The  weight  of  the   center  span  load, 
consisting  of  the  test  cylinder  and  test  cylinder  yoke,   exciter  extension 
and  coupling  hardware,    and  the  exciter  moving  element,   was  estimated  to 
be  Wc   =   5-5   lbs.      Thus,    for    0(Jt  =  2.25 

\  Wb 

—  r=-     =     5.76  from  Fig.    6 

mc  Wc 


and 


also 


W       =      5.76  x  5.5     =     31.68  lbs 


¥       =     bxhx2/x^=     31.68  lbs 


Another  equation  for  J(  in  terms  of  h  may  be  obtained  from  Equation   (l8), 

Appendix  A:  2 

k  AToO  12  V"  2T|f 

Q(         =  EIgQ  =       Egn  (    h    ) 


Then 


d  - 


12   y      (2.TTfN2 
EgQ  h 


1/k 

/       =      2.25 


78, 


Squaring  both  sides  of  this  equation  gives: 

<3±-)  (2TTf)  \-     =  (2.25) 

Ego  n 

where 

f  =  100  cps 

6    .  2 
E  =  30  x  10  lb/in 

g  =  386  in/sec 

thus  2 

,12  x  0.283, X/2  ,  x   i  N2 

( T^)     (2TT  x  100)  4—  =   (2.25) 

30  x  10  h 

.2  .2 

V  =  4T5>  or  h  =  % 

Substitution  of  this  relation  for  h  into  the  above  expression  for  W. 

b 

gives:  2 

W     =bxhx2/     x  Y   =   31.68  =  b  x  f—     x  2  /    x     V 
b  1+75 

Thus  /may  be   calculated: 

2 


.5x^x2/      x  O.283     =     31.68 


/3  ^31.68  xM5 

/  1.5  x  2  x  O.283 

/  =     26.07" 

Beam  length  =     2  Jl      =   52.  lV 
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Since  h  =  J2/kJ5: 


h  = 


(26.07) 
^75 


=  1.43' 


The  steel  beam  dimensions  for  a  resonant  frequency  of  100  cps  with  a 
center  load  of  5«5  It's  are: 

1.50"   x  IA3"   x   52.lV      . 
The   steel  beam  as  machined  had  the   following  dimensions: 
1.475"   x  1.408"    x  52.20"      . 
The   resonant   frequency  of  this  beam  was   determined  as   follows: 


Wb 

W 
c 

\ 
m„ 


O.283  x  l.W  x  1.408  x  52.20     =     30.717  lbs 
5.50  lbs 


W 


w. 


b        30.717 
5.50 


=    5.585 


From  Fig.   6,   for  m./mc  =  5.585,    0(1=  2.247;    hence: 


A     =     *       = 


/■ 


f 


2.247 
26.10 


=  0.0861 


(A 


12  Y    /2TTfy2 

Eg.  h      j 


1/4 


=      0.0861 


8o. 


12  x  0.283  ,2VTf     2  ,        0.    N^ 

(T7TO8)       =     (0.0861) 


30  x  106  x  386 

and 

f     =     97-0  cps 

It  was   desired  to  have   a  minimum  resonant   frequency  of  100  cps, 
however.      The   calculated  resonant  frequency  was   raised  to   IO5.3  cps  "by 
shortening  the  "beam  to  a  length  of   50.00".      Thus,    the   above   calculation 
procedure   is  used  again,    as: 


W  =     O.283  x  l.kjk  x  1.408  x   50.00     =     29.36T  lbs 
b 

m,  a 

—  _     29-3oT     =      5.3^0,   cfjf  =  2.2^3  from  Fig  6 

mc  5.5 

OL  -     ^^       -       0897     -         2£JC     (2TLI)2 

°  '     25.00        "      '°y9T     -  EgQ        l      h   } 


from  which 


f       =105.3  cps. 


The  actual  experimental  resonant   frequency  of  the   system  was   104  cps. 

The  nodal  distance   from  beam  end,    x   ,   was   determined  from  Fig.    7 

n' 


with  the  value  of  o(J    =     2.2^3: 

for   #/=  2.243,        ^     =     0.503 

J 
since    /    =     25. 00" 


*n        =     25  x  0.503     =     12.575" 


Each  gimbal  assembly  was   then  placed  this   distance   from  each  beam  end. 
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APPENDIX  C 
CALIBRATION  PHASE   HEAT  TRANSFER  SAMPLE  CALCULATIONS 

The   calculations   presented  below  were  those  performed  for  calibra- 
tion phase  Test  No.    9   (Table   l)  which  had  the  following  relevant  original 
data: 

Q       =45-00  watts 
m 

t       =     266. 7°F 
s 

t      =    T6.T°F 
a 


a.    t,   temperature  potential 

Formula:       At     =     t     -  t 

s  a 

t        =     266. T°F 
s 

t      =     T6.T°F 

d 

At     =     266. 7  -   T6.T     =     190. 0°F 


b.  Q  ,   power  delivered  to  test  section 
Formula:      0^     =     O.989I  Q 

Q^     =     45.00  watts 

G   =  O.989I  x  45.00  =  44.510  watts 

c.  0  ,   heat  transferred  by  radiation 

Method:      Q       =     value  read  from  Fig.    12,    corresponding  to     £  t 
A  t    =    190. 0°F 
Qy     =     0.433  watts 
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d.  Q,  ,  heat  transferred  by  free  convection  in  absence  of  vibration 

co'  J 

5/4 
Formula:   Q    =  0.01428  A  t  '      ,   watts 
co 

t  =  190. 0°F,   t5^  =  686.804 
Q    =  0.01428  x  686.804  =  9.809  watts 

CO 

e.  Qfl  +   Q,,      )    axial  heat  loss   from  test  cylinder 
Formula:      Qd  +   Q^     =     0^  -   (Qr  +   QCQ) 

q^     =     44.510  watts 

Q       =     0.433  watts 

Q         =     9.809  watts 
co 

Qd  +  \Q     =     44.510  -  (0.433  +  9.809) 
Q    +   ^      =     34.268  watts 
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APPENDIX  D 
VIBRATION  PHASE  HEAT  TRANSFER  SAMPLE  CALCULATIONS 

The  calculations  presented  "below  were  those  performed  for  vibration 
phase  Test  No.  6  (Table  2)  which  had  the  following  relevant  original  data: 

f     =  104  cps 

Q     =44.90  watts 
m 

a     =  O.O836" 

^*hm     =     52'llOF 
t_  =     228.03°F 

a 

t  =     T8.37°F 

a.   t,    temperature  potential 
Formula:  A.t     =  t_  -  t 

S      cl 

ts   =  228.03°F 

t   =  T8.3T°F 
^t  =  228.03  -  78.37  =  i49.66°F 


~b  -      Q^>   power  delivered  to  test  section 
Formula:   0^  =  O.9891  Q 
Q   =  44.90  watts 
Q   =  O.9891  x  44.90  =  44.41  watts 


Q  ,  heat  transferred  by  radiation 
r 


Method:   Q   =  value  read  from  Fig.  12,  corresponding  to  &   t 
A  t  =  l49.66°F 


Q   =  0.314  watts 
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d.  Q-j  +   0        ,    axial  heat  loss   from  test  cylinder 

Method:      Q,  +   0  =     value  read  from  Fig.    13,    corresponding  to 

^--'     hm 

A    t  =      52.11°F 

hm 

Qd  +   ^uo      =      31,63  Watts 

e .  Q       ,    heat  transferred  by  convection  in  presence  of  vibration 
Formula:      Qqv     =     ^  -   Qr  -   (%  +   Q^) 

Q       =     44. 4l  watts 
Q       =     0.314  watts 

Qd  +   ^uo      =     31,63  watts 

q   =  44.41  -  O.31  -  31.63  =  12.47  watts 

f .  h.   heat-transfer  coefficient   in  presence  of  vibration 


Formula:   h   =  3-412 


"cv 


v      "      J '  A  £  t 

Q  =     12.47  watts 

cv 

A  =     0.0993   ft2 

/\  t      =     l49.66°F 

12.47 
h  =      3.412  x  — — -7-T77     =     2.863 


v      J*      0.0993  x  149.66      '   J  hr  ft2  oF 


g.   af,  the  intensity  of  vibration 

a  x  f 
Formula:   af  =  -rr= — 

a   =  O.O836" 

f   ■=  .104  cps 

O.0836  x  104     n  „nr.   _  / 
af  =  ^— =  0.725  ft /sec 
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